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Interference Between Wing and Body 
at Supersonic Speeds—Theory and 
Numerical Application’ 


CARLO FERRARIt 


Cornell Aeronautical Laboratory 


ABSTRACT 


The supersonic flow past a body of revolution with pointed nose 
and carrying a wing is treated. The flow at infinity is uniform at 
a small angle of attack with respect to both the axis of the body 
and the chord of the wing and parallel to the plane of symmetry 
of the body-wing combination. The wing profile is assumed to be 
that of a flat plate. All perturbations on the uniform stream are 
assumed to be so small that it is legitimate to linearize the equa- 
tion of motion. 

The following are determined: 

(1) Interference of the wing on the streamlined body, assum- 
ing that the induced field generated by the wing is that which 
would exist around the wing if it were placed in the uniform 
stream alone. 

(2) Interference of the body on the wing, neglecting the distor- 
tion of the field around the body due to the presence of the wing. 

A detailed numerical example serves to show that the inter- 
ference effect of the wing upon the body is about half as great as 
the effect of the body upon the wing when the lift increments are 
compared, but the moment increment on the body is nearly two 
and a half times as large as the increment experienced by the 
wing. An optimum angle of wing setting for least drag is also 
determined. 


INTRODUCTION 


— DEVELOPMENT OF A THEORY of wing-body inter- 
ference is importunate because the effects of such 
interaction upon the aerodynamics of a supersonic con- 
figuration are apt to overshadow completely the be- 
havior that would be predicted on the basis of indi- 
vidual treatment of the body and lifting surfaces as 
though they were isolated in the impinging stream. 
Empirical data indicate that large discrepancies exist 


Received January 30, 1948. 

* Italian manuscript submitted to Cornell Aeronautical Lab- 
oratory, March 4, 1947. The numerical application was ac- 
complished under Navy BuOrd sponsorship. 
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between experiment and the simplified theory which ac- 
counts solely for the aerodynamics of the isolated com- 
ponents acting alone. Not only do the interference 
phenomena markedly alter the lift, drag, and moment 
of the combination of body and lifting surface, but, of 
course, the downwash excited by the same changes of 
the flow in the region of the wing-body juncture also 
profoundly affects the behavior of the stabilizing and 
control surfaces in general. 

The subsequent discussion is a step in the delineation 
of the mechanisms at play in the generation of inter- 
ference effects. It is based on the assumptions of a per- 
fect fluid and of the existence of a velocity potential 
(norirotational field). In addition, the method of small 
perturbations (and, consequently, small angles of at- 
tack) imposes further restrictions. The latter confin- 
ing hypothesis may be removed by use of the method of 
characteristics to solve the same problem as treated 
here. This approach has also been investigated by the 
author, and the results are expected to be made avail- 
able through the JOURNAL OF THE AERONAUTICAL 
ScIENCES at a later date. There are many other 
phases of this problem to be solved in addition to the 
analysis presented now, but the principal basic effects 
are clearly delimited for the simplified configuration 
studied. Thus these interactions are eminently worth 
consideration in order to arrive at a more lucid inter- 
pretation of the aerodynamics of supersonic aircraft. 


SYMBOLS 


Only those of most importance are included here. 


eo = potential function for the perturbed flow around 
the isolated body arising from the impingement 

of the longitudinal stream upon the body alone 

eo = potential function for the perturbed flow around 


the isolated body arising from the impingement 
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@(3) 


(x, ¥, 2) 


(x, 8) 


(é, 2, $) 


Fie. 1. 


Designation of coordinate axes, 


of the transverse component of the stream upon 
the body alone 

potential function for the perturbation field gener- 
ated around the body from the influence of the 
wing 

potential function for the combined effect of the 
above three factors in producing a perturbation 
field about the body in the presence of the wing; 
the velocity field described by this potential 


must be added to the undisturbed velocity, Vo, 
as well as to the velocities generated by the iso- 
lated wing, in order to obtain the final flow 
about the wing-body combination 

potential function of the perturbed field around 
the wing from the interference effect of the body 

coordinates relative to a system of orthogonal 
Cartesian axes oriented as shown in Fig. 1; the 
x-axis is aligned along the axis of the body and 
is directed in the same sense as the velocity of 
the undisturbed stream; the z-axis is normal to 
the x-axis and lies in the plane of symmetry and 
is directed in the sense of positive lift; the y- 
axis is normal to the plane of symmetry (along 
the wing span) and directed in conformity with 
a right-handed system; the origin is situated in 
accordance with the construction noted in Fig. 
2 for the determination of the perturbations 
about the body, but it is shifted to the wing 
leading edge in Section B 

cylindrical coordinates; the initial side of the 
angle is in the direction of the positive y-axis; 
the terminal side sweeps out a positive angle 
when it is rotated clockwise (looking in the di- 
rection of the positive x-axis) 

semispan of the wing 

wing chord 

radius of any right cross section of the streamlined 
body 

nondimensional coordinates defined as 


t=x/b, n=y/b, =8/b 


velocity of the undisturbed stream 


angle between Vo and the x-axis 
Ve = magnitude of the velocity in the undis- 


turbed stream 
angle between Vo and the wing chord 


Co = velocity of sound in the undisturbed stream 
Me = Vo/C = Mach Number of the undisturbed 


stream 


V M? — 1 = cotangent of the Mach angle 
po = corresponding value of the density in the 


undisturbed stream 


1948 


p = pj. = corresponding value of the pressure in the 
undisturbed stream 
Vi = limiting velocity (for nondimensionalizing pur- 


poses) 
= C,/Viv. = 
distance between the tip of the nose of the body 
and the parallel (circumference of a section 
normal to the body axis) at which the leading 
edge of the wing joins the body (see Fig. 7) 
x/d, n* = = 2/d 


8 
| 


> 
ll 


(A) INTERFERENCE EFFECT OF THE WING UPON THE 
Bopy 


(A-1) General Equations 


With the hypothesis of small perturbations the lin- 
earized differential equation defining the potential ® of 
the perturbed field around the streamlined body is, in 
cylindrical coordinates (x, 7, 6), 


Or? r r? 
One must determine a solution of Eq. (1), satisfying 
the condition that the derivative of ® in the direction 
of the exterior normal to the meridian line of the body 
be equal and opposite to the component along this nor- 
mal of the velocity vector for the stream flowing around 
the body-wing configuration, without interference ef- 
fects. This component is the sum of: (a) the similarly 
situated component of the undisturbed (i.e., free 
stream) velocity about the streamlined body, which 
may be set equal to V)(—_R + 8; sin @) since the slender- 
ness of the body is such as to make R of the same order 
of magnitude as 8, and where R = dR/dx; and (b) 
the similarly situated component of the perturbation 


(1) 


velocity generated by the wing, acting alone. This 
component may be expressed in the form 
Vin = Vod oF n(x/R) sin (2) 
m 


The determination of the functions F,,(x/R) will be in- 
dicated in Section A-3. 

If, therefore, one assumes ® = + + jn 
which ®‘" is the potential of the perturbation field 
corresponding to the longitudinal component of the 
stream, Vo cos 8; = Vo; is the potential of the per- 
turbation field corresponding to the transverse compon- 
ent of the stream, Vy sin 6; = and represents 
the potential of the perturbation field due to the inter- 
ference effect of the wing; then &”) and ®°) are de- 
termined, as is well known, following the methods ex- 
pounded by Karman and Moore,! Ferrari,? and Tsien.’ 


(A-2) Determination of the Potential 2 
For the determination of ®“) let 


$3 = r)r™ sin (3) 


Substituting Eq. (3) in Eq. (1), it is found that each of 
the functions ®,,“ must satisfy the equation 
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Now the solution of Eq. (4) is* 
3) = >) (5) 


where 


& = Sun x /ar f(x — ar cosh u)du (6) 


and the origin of the coordinates is at the point on the 
x-axis at which all harmonic components of the pertur- 
bation begin to arise (Fig. 2). Thus, 


= 
—a x/ar — ar cosh u) cosh u du 
r>, (3) = ( 
AB =WING CHORD a? Sercosh x/ar fo(x — ar cosh cosh? u du 
Fic. 2. Scheme of interval division for the stepwise evalua- etc. 
tions. 
The arbitrary functions f, fo, etc., must be calculated 
4 (= in such a way as to satisfy the boundary conditions, 
or? r or Ox? and therefore 
E (r@, (3) = a? — aR cosh u) cosh? u du = —F, (=) 
or = R arcosh x/aR R 
[2 @) R cosh h? ud 
fo(x — aR cosh u) cosh® u du 
etc. 
in which f is the derivative of f with respect to the param- Consider then the points P;, P2, etc., on the meridian 
eter on which it depends. line, which are the intersections of this meridian with 


The solution of the integral Eqs. (8) is performed the straight lines: x — x,’ = ar, where x,’ is the ab- 
most simply by following the procedure given by von  scissa of the right end of each of the above-mentioned 
K4érman.! The interval of integration is subdivided small intervals. The condition is imposed that Eqs. 
into smaller partial intervals, in each of which the in- (8) be satisfied at the points P; (Fig. 2). Then at the 
dividual f are considered constant. point P; with coordinates (x:, Ri), one obtains 


x1 x? x1 
aR, aR, a*®R;? ) ‘ 


where K,“) is the constant value of /; in the first small interval 


3 a® Ks (2 i+ oR? Y Ri (8’) 


where K,‘) is the constant value of f; in the first small interval 


1 3 xX) x? ( 2 xy? (=) 
4 E aR, 2 aR, 3 : R, 


with analogous significance for Ks", etc. 


Nile 


The value of K,“ can be calculated immediately In general, if x;’_1, x,’ are the abscissas of the left and 
right ends of the 7th interval, respectively, and K,“ is 
the constant value of /, in that interval, then one has, 
second; and the rest in like manner. at the point P, with coordinates (x,, R,): 


from the first of the set of Eqs. (8’); Ke“ from the 
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arcosh — — arcosh oR, oR, ( : ‘ 
= , 
Xn Xi (= =) 1] F, 
aR, aR, R,, 


etc 


From these equations the constants K,‘ can be de- 
termined successively, provided the F, are known. It 
is important to note that for the purpose of calculating 
the lift and moment of the body—not the local pres- 
sures—it is necessary and sufficient to compute only 
, and therefore only the constants K,!. 


The problem is thus reduced to the determination of 
the velocity normal to the body, induced by the wing— 
that is, to the calculation of the harmonic components 
of this velocity, according to Eq. (2). For this calcu- 
lation we let the wing be a flat plate and assume that the 
field induced by it is the same as it would be if the body 
were absent. Under this condition it is known that 
V,» reduces to 


™ — sin 6 (10) 


in the region included between the planes that pass 
through the leading and trailing edges of the wing and 
are inclined by the Mach angle with respect to Vo, 
while V,, = 0 outside of this region (Fig. 3). Conse- 
quently, 


F,=0 ifm>1 


F, = — Bo; 
for a(Ryo + Ri*) Sx SaRo+/f (11) 


where Ry and R,* have the meaning indicated in Fig. 3 
and / is the wing chord. 

For 0 < x < aR we have F,, = 0 regardless of m, 
and for aRy < x S a(Ro + Ri*) the graph defining 
the variation of V,, for a cross section is shown in 
Fig. 4. 

Letting 


Fic. 3. Regions in which the wing-induced normal velocities 
take on different formulations. 


7 
/ 
27 


Fic. 4. Wing-induced velocity normal to body in the region 
immediately aft of wing leading edge. 
aR» 
6, = arcsin ———; = — Oy 


aR ’ 
= + Op = 


one obtains, therefore, for a general function 


(12) 
2r 64 
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aR, aR,, (9) 
at 3 Xn — Xj-1 
— K;® {2 (arcosh ———— — arcosh 
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) 
n= sin @sin mo do sin sin ma do + sin 6 sin mo ao | 


then 
- sin? d@ + sin? 6d@ + ty sint 9 = 
T 
2 — 1)0 1) 
m— 1 m+ 1 
= 0, for m even 


59, — sin 20, | form = 1 > (13) 


Finally, for 


(aRo +1) S x a(Ro + R*) +1 


(see Fig. 3) the curve defining the law of variation of 


V,.» for any cross section is shown in Fig. 5. 
Putting 


§, = arcsin (x — aRy — 1)/aR; On = — On 
Ber = + Op = 2r 64" (14) 


one obtains 


Fic. 5. Wing-induced velocity normal to body in the region 
immediately aft of wing trailing edge. 


ala 


Fy 


Il 


3 
Il 


sin (m + 
m+ 1 


sin (m — _ 
m— 1 


0, if m is even 


— Be 6D’ — Ba 
6 d@ + sin? 6 d@ | = ——| — 264 + sin 204° |, form = 1 
6a’ 6c’ 
— Bo 
| sin 6 sin m@ dé + sin sin mé ao] 
6a’ 6c’ 


(15) 


| if m is odd and > 1 


In the event that the body diameter is so large that 
| < aR,*, then Fig. 3 no longer applies, and a new form- 
ulation for the V,, is required in the interval 


aRo +1 <x < + Ri*) 


The above-given equations must be replaced by 


F, (sin 264 — sin 204) | 
form = 1 
(m — — sin (m — 1)0y _ 
m—1 
sin (m + 1)@, — sin (m + ad 
m+1 


for m odd and > 1 
= 0, for m even 
and Eqs. (13) and (15) apply to the regions 
aRo aR +! 


and 


a(Ro + Ri*) <x 
respectively. 


S a(Ro + +1 


(A-4) Determination of the Influence of the Wing on the 

Lift and Moment of the Body 

In order to calculate the overall interference effect 
of the wing upon the body—i.e., the variation of the 
resultant of the pressures acting on the surface of the 
latter, not the local pressures—the evaluation may pro- 
ceed from the approximate solution of the first equation 
of Eqs. (8). The approximation is obtained for the 
limiting case wherein the radius of the body tends to 
zero; this evaluation is valid because only f; is of im- 
portance as far as the overall lift and moment are con- 
cerned, as pointed out above. 

Now, if we let x’ = x — aR cosh u, the middle ex- 
pression in the first equation of Eqs. (8) reduces to 


IR 


= x—aR (x x’)? 
fle’) V (x — x’)? — 


‘ 
: 
aN, 
Hi 
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if one takes the limiting value of the integral as R > 
0.’ Setting this expression equal to the right-hand side 
of the first equation of Eqs. (8) and differentiating 
with respect to x, one obtains 


fi(x) = (d/dx)[R?Fi(x/R)] (16) 


On the other hand, the variation of lift of an element dx 
of the body, produced by the velocities normal to it 
that are generated by the wing, is given by 


d RF, “*) 
= (Aula) over sin = 
0 
= | fi(x — aR cosh u) X 
arcosh x/aR 
cosh u au 


pVorrfi(x) 


again with the hypothesis that R tends to zero. The 
overall change of lift due to the field of normal velocities 
is then, taking into account Eq. (16): 


Ail, = filx) dx = 
0 


x 
2 — 2 — = 
pVo rf F, ( dx 


because F\(x/R) is continuous in the interval and zero 
at its ends. Moreover, since the change of lift due to 
the velocity components along x, induced on the body 
by the wing, is given by 


2p Vo? 
Ad, = RB.(x, 0) sin 6 d6 dx 
a 0 0 


and, by making use of the relations already given in 
Eqs. (11), (12), and (14) for determining when (x, 6) 
is zero or equal to a constant, 2, one perceives that 


a(Ro+ Ri*) x — aR,\2 aRot+l 
tale = - (= 2") Rdx 
Qa aRo aR a(Ro+ Ri*) 


Even if R is finite but constant, this expression reduces 
to 


= = ALy (17) 


so that this change is equal to the lift of a portion of 
wing which has just the same extent as the wing area 
covered by the body. When R is extremely small, it is 
necessary to observe simply that 


d(AgL,)/dx = (4/a)pVo2BoR 


for x lying between the abscissa x» = aRp (of the leading 
edge) and the abscissa x) + / (of the trailing edge), 
while d( A,L,)/dx = 0 for x outside the interval. Con- 
sequently, the overall change of lift of the body is 


AL, = AcLy = Rdx 


With the same approximation we may obtain the 
changes in the moment, about the y-axis, of the aero- 
dynamic forces acting on the body as a result of the 
interference effect of the wing. This change is given 


by 
0 1 R 


Under the hypothesis, generally quite legitimate, that 
R may be considered constant for that portion of the 
body influenced by the wing, one gets 


= — pV (18) 


a(Ro+ Ro*¥)+1 2 
aRot+l aR 


The amount of moment change, A:,M/,, due to the 
longitudinal velocity induced directly by the wing, may 
again be obtained simply in terms of a finite R, provided 
it is constant in the region 


aRy + R*) +1 


Proceeding as in the case for Eq. (17), it is found that 


4 1 R 
= -(*) Bo () + 


and this moment is likewise, when R becomes very 
small, exactly equal to the moment that would arise 
from an isolated portion of wing, the spanwise extent of 
which is 2R. 


(B) INTERFERENCE EFFECT OF THE Bopy ON THE 
WING 


(B-1) Condition That Interference Be Nil 


The perturbation potential field for an isolated body, 
considered immersed in a uniform stream of velocity 


Vo, is 
4+ 2) = Vol x/ar f(x — ar cosh u) du — 


a sin x/ar — arcosh u) cosh udu} (19) 


At any point of the (x, y)-plane, the velocity, corre- 
sponding to this potential, has a component normal to 
that plane given by 
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(x — ay cosh u) cosh udu (20) 


0 
Qa 
arcosh x/ay 


wherein one may set, at least in first approximation,’ 
g(x) = 


where S is the area of a general cross section of the 
body. 

Within the degree of approximation consistent with 
the linearized theory, the condition for the nonexist- 
ence of body interference effects on the wing, when the 
latter is at a given angle of attack, 62, is that the 
slope of the wing profile at each point in the plane 
of the wing [considered to be close to the (x, y)- 
plane] be 


dz 
g(x — ay cosh u) cosh u du (21) 
dx Y Sarcosh x/ay 


Under this restriction, certainly, the velocity induced 
by the body has a zero component normal to the wing 
profile. 


(B-2) Condition of Non-Zero Interference—General Equa- 
tions 
If points along the chord line of the profile do not 
satisfy condition (21), the velocity component normal 
to the wing [that is, to the (x, y)-plane, since the wing 
is considered to be close to the (x, y)-plane], at each 
point of its surface, may be expressed in the form 


> Gal8) cos (22 


Here the origin of coordinates is at the leading edge of 
the profile which lies in the (x, 2)-plane, and the func- 
tions G,,(~) can be determined from Eq. (20) as indi- 
cated in Section B-5. In order to satisfy the boundary 
condition relative to the wing, it is necessary to 
superimpose on potential field (19) a new set of ve- 
locities with a corresponding potential that is ex- 
pressed as 


mn’ cos (xmn/2) (23) 


and which is defined by the following conditions: 
(a) it must satisfy 

(0°b’/On?) + = (24) 
everywhere; (b) while in the plane of the wing with 


¢ = 0, it is required that 
Gn(§) (25) 


By substitution of the potential (23) into Eq. (24), one 
obtains 


cos + cos 


where k = rm/2; so that it is apparent any ®,,’ must 
satisfy the equation 


while Eqs. (23) and (25) give, on the plane of the wing 
=0, 


= —Gn(E) (27) 


Now Eq. (26), with condition (27), is formally analo- 
gous to the so-called telegraph equation’® (in its reduced 
form) and, in consequence, can be solved by means of 
the corresponding Riemann functions. In order to 
give the procedure a more well-rounded development 
than was done under Section A, however, the solution 
of Eq. (26) that satisfies condition (27) can be found as 
follows: 


(B-3) Determination of the Potential ®,,' 


(B-3-1) Construction of the Potential from Elemental 
Sources, and Related Doublet Distribution—We now 
consider ®,,’ as being a function of » also, and satisfying 
the equation 


A particular integral of Eq. (28) is immediately seen to 
be 


= [cos (k/a)r]/r (29) 


in which now r = Ve — + £°). The function 
®,,’, as defined by Eq. (29), is called the potential of the 
source that is located at the origin of coordinates. If 
the source, instead of being at the origin, were located 
at the point with coordinates, 


the corresponding potential would still be expressed by 
Eq. (29), in which now, however, 


Consider further a uniform distribution of sources on 
the semiplane (£, > 0) extending from 7’ = 
—« ton’ = +; the intensity of source in the ele- 
ment dt’ dn’ is h(£)dé’dn’. The potential of the source 
distribution will still satisfy Eq. (28), since it is a linear 
combination of particular solutions of this equation. 
Consequently, we have: 
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h(t’) cos (k/a) VW (E — — a2[(n — 9’)? + dt’ 


cos (k/a) V(E — £')? — a2 [(n — 0’)? + dn’ 


/a) 


(30) 


VE — — a [(n — 9’)? + 


Making the transformation 


= (1/a)V (é — §’)? — at cos ¢ 


and since 


—dn’! = —(1/a)V — £’)? — sia 


VE — &)? — — + = 
— 


— a? sin g 


one arrives at 


1 
X 
a Jo 
k 
[ cos VE af? sin e| dy 
0 a 


|! at | ay 
(31 


where J) is the cylindrical Bessel function of zero order, 
This result satisfies Eq. (28), and therefore, since it 
does not contain 7, it also satisfies Eq. (26). 

Corresponding to this potential for the source dis- 
tribution there will exist a potential for a doublet dis- 
tribution,’ with axis in the direction of ¢ and lying in 
the plane (£, 7), of the form: 


a OF Jo 


k 
= + rh[t — al¢|] Ji [Eve | 


where J; is the cylindrical Bessel function of the first 
order and where the minus or plus sign should be used 
according to whether ¢ is greater or less than zero, re- 
spectively. It is seen from Eq. (32) that the deriva- 
tive of ®,,’ with respect to ¢ is continuous at the plane 
¢ = 0, while that with respect to £ is discontinuous and 
has the value 


(B-3-2) 
as to satisfy Eq. (27). 
the result that 


0 
= tahn(~ — at) — hf 


— cosh u) cosh u J; (k sinh u) du — 
arcosh 


| k 


dé’ 
VE = = 


(32) 


In what follows we shall consider the perturbed flow in 
the semiplane ¢ > 0; the remark just made above 
enables one to determine the flow in the other semiplane 
when that in the first is known. It is assumed that 
®,,’, defined by Eq. (23), is to be expressed by Eq. (32), 
into which it transcribes when /(é’) is replaced by 


hn(&"). 


A pplication of the Boundary Condition.—The function h,,(’) a priori arbitrary, must be determined so 
In order to calculate 0®,,’/O0¢, it is convenient, as usual, to put t’ = § — af cosh u, with 


hy(& — af cosh u)Ji(Rk sinh u) du + 


0 
hy(& — af cosh u) X 
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= — at) + hy Ji [Eve - | — — a%? 


k? 
hnlt!)Si [* ~ ~ | dé’ 


where h,, and J; are, respectively, the derivatives of h,, and of J; with respect to their corresponding parameters. 
On the plane of the wing, for which ¢ = 0, one has 


a Jo a a 


0 Qa a Jo a 
but J:(6)/@ + J,(@) = Jo(@) and therefore 


+ hal’) Jo| | ae 
rk f dt’ = (33) 


(B-3-3) Stepwise Procedure for Evaluation of the Doublet Distribution.—Equation (33) can be solved simply by 
following a procedure analogous to that followed in Section B-2. If the interval 0 to //b (where / is the length of 
the wing chord) is divided into parts, in each of which h,,(£’) is assumed to vary linearly, then 


= 


in the first small interval, 
= Ai + Hy — As) 
in the second small interval, 


hn(!) = di + + Hm — — 


in the third, and so on. 
If, for the sake of simplicity, the lengths of the intervals are all set equal to A, then one has in the ith small in- 


terval that 
= + Hele’ — @ — 1) 4] 


ref 


k k 


and 
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Therefore, imposing the condition that Eq. (33) be satisfied at the point § = A, cemeapeneng to the right-hand 
end of the first small interval, one obtains: 


H, | (: a) + a) + 2nkA (*)) = —Gn(A) (34) 
a 


whence /7,,“) can be determined. 
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Analogously, imposing the condition that Eq. (33) be satisfied at the point £ = 2A, corresponding to the right- 
hand end of the second small interval, and since 


§ Jo E (: )| dt! = + Jy dt! + f° — A) Jo dt? 
0 0 


while 
2 
0 Rk? la a a a 
2 —A k 
and 
@ 
J Jy dé’ = > Juss] (: 
A n=0 
and 


and: lastly 
dt’ = So* Ji dt! + Hn® Ji dé’ 


then there evolves, for & = 2A: 


kA kA A 
Har (* =*) + ) +2 = —G,,(24) 
or, taking into account Eq. (34), 


Hye (= *) + 4kA Juss + aJo 
H,, x () + kA E + 2 (**)]! = —G,,(A) — G,,(24) (35) 


for which H,,“ is determined. Similarly one obtains 


BRA 3kA BRA 
EX (=) + 6RA Juss (7 ) + Ga + + 
a a 


2kA 
4kA Juss *) + aJo (= + H,, [eas + 2kRA DY (‘) + (**)] = 
Qa a 


—Gn(A) — Gn(24) — Gn(3d) (36) 


which gives //,,“, and so forth. 


(B-4) Solution of the Problem in the Limiting Case of an Infinitesimally Narrow Chord Wing 
An approximate solution of Eq. (33) can be obtained by taking as the solution of that equation the one corre- 


sponding to the case of an infinitesimally narrow chord wing. Indeed, substituting in Eq. (33) for Jo and J, their 
limit values, when their argument tends to zero, Eq. (33) reduces to 


rathm(£) + 1(k?/2a) So hm (€’)dé’ (33’) 


In order to derive this Eq. (33’) from Eq. (33) and for the purpose of demonstrating within what limits the sub- 
stitution of Eq. (33’) for Eq. (33) is legitimate, the following detailed procedure is subjoined: 
The familiar expressions for the Bessel functions are: 


Toth — &')/a]} = 1 — — + . 


and 
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and 


These series are uniformly convergent in every finite interval; hence, the integral of each series in Eq. (33) may be 
replaced by the following series of integrals: 


Jo E dé’ hn(é’) dt’ — aL inte’) di’ +. 


Now let = x/b = (x/J)(1/b) = (1/b)X and &’ = (1/b)X’; then 


0 Qa 0 0 . 


where O(k*/?/4b?a*) is an infinitesimal of the same order as k*/?/4b?a?. Analogously, 


han 8’) E (= dt’ = ) thes 
(X’)(X — dX’ + 


and 


Therefore we have, 


If one assumes, accordingly, that k//2da is an infinitesimal of the first order, it is legitimate to neglect the term 
O(R*/?/4b?a*), and Eq. (33) becomes Eq. (33’). 

It follows, because of the above discussion, that in the case of a given wing chord for which //b < 1 the approxi- 
mation supplied by means of Eq. (33’) improves as a becomes larger and as k (that is, m) becomes smaller. For 
large values of m the approximating Eq. (33’) can still give a close estimate of the value of h,,(£) only at the points 
P of the wing which are in proximity to the leading edge, O; they must be so situated that, if ] is set equal to OP, 
then 


= <1 
Now consider the approximating Eq. (33’) further. By differentiating it with respect to £, we obtain 


+ = —Gn(E) (37) 
while at the same time it must be true that 
h,(0) = 0; hm (0) — (1/ra)G,,(0) (38) 


One deduces that the solution of Eq. (37), which satisfies the conditions (38), is given by 


= cos (5) f° Gn(é) sin (5) dé — 


sin G,,(£) cos dt — G,,(0) sin (39) 
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Fic. 6. Antisymmetrical extension of body-induced velocity 
distribution. 


(B-5) Determination of the Harmonic Components of the 

Perturbed Velocity Induced by the Body 

The solution of the problem is thus reduced to the 
determination of the functions G,,(£) which represent 
the amplitudes of the harmonic components of the 
velocity, perturbed by the body, at the wing. The 
calculation of these functions, when the potential of the 
perturbation field induced by the body is known, can 
be performed in the following manner. Once the graph 
of the normal perturbation velocities V,,,,/Vo has been 
determined along the span of the wing for a given value 
of £, by means of Eq. (20), it will be possible in general 
to approximate the function defining V,,,./V» with as 
parameter (for » varying in the interval from —1 to 
+1) by a trigonometric polynomial of period 4 and 
consisting of terms. 

If the graph of the function V.,,,/ Vo is extended anti- 
symmetrically (Fig. 6) outside the interval -1 <7 < 
1, only the terms that contain odd values of m will be 
different from zero (that is, m = 2n). 

If (Vaw/ (Vnw/ Vo)i, are the values of Vo, 
respectively, at the points 


(7 = 0; & = &) 
and 
E ~ for (i 
s4 ¢ then the value of G,,(£,) can be calculated from the 
formula of Clairaut de la Vallée Poussin :’ 
min 
Gn, =- — 
V, 
nw —1)" nw 
( ae) + ‘ 
2n 
form = 1,3, ...,2" — 1 (40) 


. oy (B-6) Calculation of the Changes of Lift and Drag of the 
er Wing Caused by Interference Effects of the Body 


The component along £ of the perturbation velocity 
corresponding to the interference effect of the body on 
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the wing is, at the plane of the wing, the sum of: 
(a) a component corresponding to the potential 6“) 
and given by 


u, = Vo {an x/ay fx — ay cosh u) du 


where the abscissa x must be measured from the apex 
of the nose of the body; (b) a component correspond- 
ing to the potential 6’ just determined; it is given by 


us = 06'/Ox = Vo cos (xmn/2) 
Of 
= —Vor x cos (1mn/2)hm(E) 


for the upper surface of the wing, and has the same 
value but of opposite sign for the lower surface of the 
wing. 

It follows that the pressure change on an element 
b*dé dn of the wing surface is equal to Apg = —pVy X 
(uw; + ue) for the upper surface and Ap, = —pV» X 
(u; — U2) for the lower surface, where p is the density 
of the fluid. Thus the change of lift on such an element 
which is caused by interference from the body is: 


@AL = b*déidn( Ap, Apa) = 2p Voueb*dé dn 
= cos hin(€) dé dn 


At the wing tip, for » = + 1, and m being odd, one 
always has d*AL = 0. Moreover, 


dAL (;) 


The overall change of lift is therefore 


AL = cos (wmn/2) dn 


8p Vo? (- 1) (41) 


In order to keep the wing lift unchanged, it is thus 
necessary to change the wing incidence an amount 
ABs given by 
ALa AL 
~ 4blp Vo? 

Hence, the increment of wave drag (induced drag), for 
the same lift, is 


AL 
= 


kes 


AB: 


AD 


where L,, is the wing lift. Therefore, the increment of 


drag coefficient is given by 


ACp = 2a (—1)*? [hn(l/b)]/m (42) 


in which C;, is the lift coefficient of the wing. 
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(C) NUMERICAL APPLICATION TO DETERMINE THE 
INTERFERENCE OF WING UPON THE BoDy 


(C-0) Description of the Wing-Body Configuration and 
the Conditions of the Free Stream 

A numerical application has been carried through for 
the case of a wing-body configuration having the fol- 
lowing properties: 

Assume the origin of coordinates to be at the tip of 
the nose of the body, and let \ be the distance between 
the tip of the nose of the body and the diametral cross 
section at which the leading edge of the wing joins the 
body (see Fig. 7). In terms of the nondimensional 
coordinates &* = x/d, n* = y/d, and = the 
body contour is delimited by means of the relations: 


nt= — in the interval 0 < <1 (43) 


n* = 0.1 = Ro/d, for 1 < &* < 1.3 (43’) 
n* = 0.068 + 0.049&* — 0.01886¢*?, 
for 1.3 < &* < 1.75 (43”) 


The wing is rectangular and is located in the plane 
= 0 (see Figs. 7 and 8); its chord lengthis/ = 1.4Rp, 
and the semispan is b = 4Ry = 2.86. 


— 


Fic. 7. Wing-body configuration for which the interference 
effects are calculated by numerical application of the small per- 
turbation theory. 


The conditions of the undisturbed flow under con- 
sideration are defined as follows: 

The angle of attack of the body with respect to the 
x-axis is 8; = 8°. 

The Mach Number is /, = 2, with the result that 
a= — 1 = 1.732. 

(C-1-0) Calculation of the Perturbation Potential 
About the Body.—In order to counteract the potential 
field related to the velocities induced by the wing, 
which are normal to the surface of the body, the distri- 
bution of singularities along the body axis which ac- 
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Distribution along body of increment in lift produced by velocities normal to the body, induced by the wing. 
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complishes this is now determined. Take the origin of 
-axes to be located at the point O;, which is the beginning 
of the portion of the axis throughout the extent of 
which the distribution is built up (see Fig. 8). It is as- 
sumed that the supersonic doublet strength is constant 
over each segment of the x-axis of length A* = A/d = 
0.05. 

(C-1-1) Calculation of the Constants K®.—From the 
first of Eqs. (8’) one obtains the following, in the case of 
point P; (located on the meridian line that belongs to 
the plane 6 = 7/2), which possesses the coordinates 
&* = 0.2232, ¢* = 0.10: 


2: 2.23\? 
; (3) Ky (arcosh ( =) 1) = 


1.73 
—2 in 20 
2 


Now we have 
sin 0, = — = 0.289 
and, consequently, 
F, (2.23) = —0.00146 


From these values we find that A,“ = —0.000545. 

Analogously, it is found that the next value of K,‘” 
is K,; = —0.00442. 

In determination of K,‘ it should be borne in mind 
that neither Eqs. (13) nor (15) are applicable, but it is 
necessary to employ the expressions involving both @, 
and 6,4 given at the end of Section A-3. With this 


caution it is found that K,“. = —0.01709. 

In similar manner, the successive values of the K,‘” 
are computed to be: AK, = —0.03652; K,® = 
—0.00631; K,® = +0.02124; AK,” = +0.03636; 
K,® = +0.01701; K,® = +0.00868; K,0 = 
+0.00399; = +0.00112; = —0.00068; 


= —0.00177;, = —0.00238; K," (at 
the end of the body) = —0.00266. 

The calculation of the A,“ values for the points 
located on the cylindrical portion of the body goes for- 
ward rapidly because many of the constants used in the 
previous steps may be re-employed. In the case of the 
points located more to the rear on the body, however, 
the section diameter progressively decreases, and con- 
sequeritly all the terms must be recomputed anew at 
each step. 

The resulting pattern of variation of the K,‘? values 
along the x-axis is depicted in Fig. 9. 

The values of K,,“”, for n different from unity, have 
not been determined because they are of no interest 
in the calculation of the resultant overall action at each 
diametral section (parallel) of the body. 


(C-2) Evaluation of the Change in Lift on the Body Which 
Is Due to the Normal Velocities Induced by the Wing 


The field of normal velocities, induced by the wing, 
produces a change in lift on an element of length, dx, of 


1948 


the body. This increment (or decrement) is given by 
the relation: 


[ 1 
Ven? ILd(x/Ro) Ry 


0 
I Ji (x — @R cosh u) cosh u dy 
arcosh x/aR 


From this formula we now write: 


B 1 ia 


which applies to the parallel with abscissa x, and radius 
R,. Thus we have: 


For the parallel with abscissa &* = 0.223: 


‘B, = 2(1.732)(—0.000545) V (2.23/1.73)? — 1 = 


— 0.0015 
For the successive parallels we find: B, = —0.0132; 
B; = —0.0550; By = —0.1330; Bs = —0.0960; 
Be = —0.0193; Br; = +0.0584; Bs = +0.0540; 


By = +0.0480; Bry +0.0413; Bu = +0.0346; 
By = +0.0282; Bi; +0.0224; Bu = +0.01715; 
By = +0.0126. The calculations in this case are 
again easily accomplished because of the reuse of pre- 
viously calculated constants from step to step; how- 
ever, after Ps with abscissa = 1.3 based on the system 
of coordinates with origin at the nose, the progressively 
decreasing diameter of cross section requires the re- 
computation of each term in the equations at each point. 

From these values of B; there evolves the distribution 
of lift along the axis of the body which is depicted in 
Fig. 8. The change in lift A,Z, may be calculated by 
means of the summation: 
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Fic. 9. Stepwise variation along the body of the supersonic 
doublet strength. 
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INTERFERENCE BETWEEN WING AND BODY 


— 0.00074 


If the wing surface, 20/, is taken as the reference area, 
the lift coefficient that corresponds to AL, turns out 
to be 


= = —0.000207 


(C-3) Computation of the Change in Aerodynamic Mo- 
ment of the Body Caused by the Normal Velocities In- 
duced by the Wing 
The change in aerodynamic moment of the body, 

taken about the leading edge of the wing, which is due 

to the normal velocities induced by the wing, is defined 


as 


AiMyp = x [d( AiL,)/dx] dx 


where the integral is taken over the whole extent, L’, 
of the body axis for which values of d (A,L,)/dx dif- 
ferent from zero exist. Since the values of d( A,L,) /dx 
may be assumed to vary linearly between the consecu- 
tive parallels of the body, the above indicated integra- 
tion may be written as 

1 


1 
My = -- — 


Xi + Xi B; 
( ) | 


The result of this calculation in the case under consider- 
ation is found to be: 


x 


V,,2] 41 My = —0.4268 


If the reference volume is taken as 20// instead of 7R,’, 
then we obtain 


—0.4268 ( )- — 0.0855 


1 


(C-4) Determination of the Change in Lift and Moment 
on the Body Which Is Due to the Longitudinal Compo- 
nent of the Induced Velocities 


Since the diametral cross section is constant as far 
rearward as Ps (with abscissa = 1.3), because a(Ry + 
R,*) + | reaches back only slightly farther to station 
1.313\, and because the body radius only reduces to 
R = 0.998R, at the next calculated point (P7), it is per- 
missible to use the following simple formula for the 
increment in lift on the body: 


(D) 
(D-0-1) 


—In order to determine this velocity component one needs to calculate the function g(x). 
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Fic. 10. Body-induced velocities normal to wing. 


4 
(V,,2)(0.1396) (2R) (2) 


= (2.31)(0.1396)("/s) = 0.0806 


based on the area 2)/. The total change in lift on the 
body due to interference from the wing is, thus, 


Cy, = —0.0002 + 0.0806 = +0.0804 


For the case of the moment change due to the longitu- 
dinal component of the induced velocity, we use the 
formula: 


= (6) 2R) (- ) 
a 2 4+ 
Thus, the change in moment coefficient, based on the 
volume 26/?, is: 


4 R R 
= — + = (a) (*)| 


—0.322(1/s)(1 + 1.943) 


—0.119 


When this increment in pitching moment is added to 
the moment due to the normal component, the total 
change is found to be: 


ACup = —0.1186 — 0.0855 = —0.289 


NUMERICAL APPLICATION TO SHOW THE INTERFERENCE OF BODY UPON THE WING 


Calculation of the Component of Velocity, Induced by the Body, Which Is Normal to the Plane of the Wing. 


This function is pro- 


portional to the unit intensity of the doublet distribution located along the axis of the body (which is assumed iso- 
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lated at present), and it is referred to the transverse velocity, V.. sin §;. , If S denotes the area of an arbitrary 
general parallel, then this function is given approximately by the formula 


g(x) = (B:/)(dS/dx) 
Taking Eq. (43) into account, one obtains 


g(x — ay cosh u) = '/58,A[2&* — + + (—2 + BE* — 3E**)an* cosh u + 
(—3 + 3&*)a’n*? cosh? u — a*n** cosh® uj (44) 


Now since &* — an* < 1 for all the points of the wing under consideration, it follows that the component of veloe- 
ity, induced by the body, which is normal to the plane of the wing is given everywhere by 


0 
= g(x — ay cosh u) cosh u du 
atcosh £*/an* 


an* an* a’n 
sinh + arcosh 2 an* Vay 1+ sinh 4u* (45) 
wherein we have set 
A, = 2&* — 4+ A, = —2+ A; = —3 + 3&* (45’) 
while 
sinh 3u* = — 1 — 1] 
and 


sinh = V (E*2/a2n*?) 4(*/an*) V (E*2/02n*?) 


For the selected values of 6; and a in this problem, the following values of v.o/v,, are obtained at points on the wing, 
by use of Eq. (45): 


At the Leading Edge: At the Wing Axis: At the Trailing Edge: 
= 1 &* = 1.07 = 1.14 
n* 0.1 0.2 0.3 0.1 0.2 0.3 0.1 0.2 0.3 
2 0.151 0.0400 0.0164 0.152 0.0417 0.0182 0.150 0.0426 0.0196 


These values of v/v, are plotted in Fig. 10. Since the differences between the values of v/v, at each n*, for 
any of the three values of &* considered, are extremely small, a single diagram has been selected as defining the 
law of variation of the velocity components normal to the plane of the wing. This diagram is shown in solid line 
on Fig. 10. 

(D-0-2) Determination of the Harmonic Components of the Velocity Perturbation Produced by the Body.—Through 
use of the graph of v.9/v,, presented in Fig. 10, the harmonic components of the velocity; perturbed by the body, 
which is directed normal to the plane of the wing are now computed. The formula of Clairaut de la Vallée 
Poussin may be written in this case as: 


7 m1. 
Gm => : -> (*) cos + *) (—1) ls (for 1, 3, 5, 7) (46) 


where we have let (v.0/v,,)y) and (v,0/v,,); be the values of at the points 7 = 0, and n = 0.257 fori = 1 to8 
These coefficients describe a function v/v, which is derived from the graph of Fig. 10 by extending it antisym- 
metrically with respect to the straight line » = 1 and which is then approximated by a trigonometric polynomial 
of four cosine terms. 

Since we have: 


n= 0 0.25 0.5 0.75 1 1.25 1.5 1.75 2 
= = 0.151 0.151 0.041 0.018 0 —0.018 —0.041 —0.151 —0.151 


it follows that G, = +0. 1254, Gs; = +0.0438, G; = +0.0027, and G; = —0.0210 The curve defined by Eq. (22), 
with G,, values just derived, is depicted in Fig. 10 by means of a dashed line. 
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D-1) Calculation of the Potential for the Perturbation About the Wing 


In order to counteract the component of velocity, induced by the body, which is normal to the wing, a super- 
imposed velocity field is to be determined. This is accomplished by dividing the chord of the wing into three 
small intervals of equal magnitude, A. In each of these intervals it is assumed that the function h,,(£), contained 
in the expression for the potential of the velocity field, will vary linearly. 

(D-1-1) Calculation of the Constants H,“.—The constant value of dh,/dé in the first small interval is denoted 
by i‘. Now in this case A = //3b = 0.117; k = 7/2; and therefore (k/a) A = 0.106; it follows that: 


Ee 73)Jo(0.106) + (5 ) (0.117) J;(0.106) + (5 Jo. 117 Jon41 (0. 106) | = —0.125 
But we have:? Jo(0.106) = 0.997; J,(0.106) = 0.053; while all of the Jo, + (0.106) are negligible, so that JZ,‘ = 


—0.0230. 
In an analogous manner, we proceed by letting /7,°) denote the constant value of dh,/dé in the second small in- 


terval, and we find that /7,°) = —0.0228, by making use of the fact that G;(A) = G,(2A) in the case under con- 
sideration, because G,, is independent of £ here. 
Likewise we find that /7,“ = —0.0225. In the latter two calculations it is also found that the Jo, + 1(0.212) 


and Jo, + 1(0.317) are negligible for m = 1 or more. 

(D-1-2) Calculation of the Constants H;“.—Now in this case k = 37/2; RA/a = 0.317; and we find that 
H; = —0.0079. 

For the value of 773“ it is necessary to evaluate J3(0.635), which may be accomplished by use of the recurrence 


formula 


+ Jn 1(x) = (2n/x) Jn (x) 


Here we have J)(0.635) = 0.902 and J,(0.635) = 0.302, and thus J2(0.635) = 0.049, and then J3(0.635) = 0.005. 
With these values of /,, we find that /73°°) = —0.0073. 

In the case of H;‘* it is necessary to evaluate J3(0.952), but the other J2, + , form > 1 are disregarded. It is 
found that 73°? = —0.0064. 


(D-1-3) Calculation of the Constants H,‘.—Here k = 52/2; kRA/a = 0.529, and we find that 17; = —0.00046; 
= —0.00037; and = —0.00025. 
(D-1-4) Calculation of the Constants H;‘’.—Here k = 77/2; RA/a = 0.741, and we find that H7;" = +0.00339; 


H;? = +0.00223; and H;“ = +0.00085. 

(D-1-5) Approximate Determination of HH,“ by Means of Eq. (39).—It seems to be in order to compare the 
results just obtained with those that would be computed by means of Eq. (39) in the case of an infinitesimally 
narrow chord wing. 

Since G,, is independent of £ in this example, Eq. (39) reduces to: 


h,,(&) —(Vf2 ‘kr)G,,(0) sin (kE/av/2) 
and thus we know that 
dh,,(&)/d& = (1, ra)G,(0) cos (RE/av/2) 


Then if we assume that the values of //,,“” are those of dh,,/d& which obtain at the midpoint of each of the three 
small intervals previously considered, it is found that: 


= —0.0230; = —0.0229; = —0.0227 
= —0.0080; = —0.0076; = —0.0068 
= —0.00048; = —0.00042; = —0.00029 
= +0.00372; = +0.00272; H;® = +0.00100 


Itappears from the above that the agreement between the results obtained by recourse to the approximate formula 
and those computations previously carried out is fairly good for the smaller values of £. In fact, in many cases 
the approximate solution should be sufficiently accurate for ordinary purposes, with a consequent great saving in 
time. 


(D-2-0) Lift Distribution Along the Chord for the Root wise lift distribution for the center profile (7 = 0). 
Section —The above-determined values of H,,“ make Let Ap denote the pressure difference between the 
it possible to calculate without further ado the chord- lower and upper surface of an arbitrary general profile 
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Consequently, it is found that: 
though the motion were two-dimensional, comes out 
Ap/('/2)p,.. = 0.3507 (in the first small interval) tates 


= 0.3546 (in the second small interval) 


= 0.3551 (in the third small interval) Ap ( 4 \2) 
= == — ) = 0.34 


When these data are plotted, the graph of Fig. 11 
results. It is worth while to note that the value of This close figure is not surprising because the value of 
Ap/(*/2)p,. V..2, calculated by use of the value v/v, v/v, is almost constant over a fairly wide range of 7 
which exists at the center section (root chord), and as_ (see Fig. 10). 


(D-2-1) Spanwise Lift Distribution.—The spanwise lift distribution is defined by: 


Now since 47(b/l) = 35.9; while 4,(//b) = —0.00796, h3(//b) = —0.00251, h;(//b) = —0.00013, and h,(//b) = 
+0.00076; and, if we let F(n) = >> h,,(l/b) cos (xmn/2) and let F\(n) = —42(b/l)F(n), then the following table 


may be constructed: 


n ad 0 0.1 0.25 0.3 0.4 0.5 0.75 
F(n) = — 0.00984 —0. 00985 —0.00897 —0.00814 — 0.00577 —0.00323 —0. 00118 
Fi(n) = 0.354 0.354 0.322 0.292 0.207 0.116 0.0406 


These data determine the graph of spanwise lift distribution which is depicted in Fig. 12. It is appropriate to 
compare this graph with the one that would be obtained if the lift at each profile section along the span were cal- 
culated by means of the equation for the two-dimensional case and based on the use of the value of the local angle 
of attack, v/v... The ordinate values of this curve would therefore be given by 


Fi*(n) = (4/V M2 — 1)(0.0/v4) = 2.309 (v.0/v,) 


(in the example being investigated). This second curve, which is a proportionally scaled-up representation of 
v~/V,., iS also drawn in on Fig. 12. The difference is, on the whole, not appreciable. 
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INTERFERENCE BETWEEN WING AND BODY 


(D-2-2) Computation of the Change in Total Wing Lift.—The total change in lift of the wing is given by Eq. 
(41) and is now evaluated to be: 


Thus the total change in lift coefficient of the wing-body configuration is: 
C, = 0.166 + 0.080 = +0.246 


This last calculation shows that the increment of lift on the body is about half as large as that on the wing (the 
latter has included an induced effect in the region occupied by the body). Actually, however, the increase of 
lift is probably less than that just calculated because: first, the lift of the body produced by the longitudinal 
velocity component induced by the wing is smaller than the lift on that portion of the wing which is covered by 
the body; second, the stream impinging upon the,wing, considered as a whole or sectionally, has a velocity less 
than V,, as a result of the longitudinal field of flow about the body. This alteration of the flow is not taken into 
account in the linear theory. 


(D-3) Calculation for the Change in Aerodynamic Moment of the Wing, Taken About the Leading Edge 


The aerodynamic moment of the wing which is induced on it by the presence of the body is (where the leading 
edge of the wing is assumed as the axis about which the moments are taken) : 


+1 
m -1 0 
(m+ 1)/2 A? 
= — }(-1) \>)— > + 8Hn® + Hn) 

m mr b b 3 
and upon substitution of the previously obtained values of H,,“? it is found that AM,/(!/2)p, V..2(b*) = —0.0208; 
and, if the reference volume is taken as 20/?, we have ACyy». = —0.0848. Then the total change in the moment 


coefficient, taken with respect to the leading edge of the wing, is ACy, = —0.085 — 0.204 = —0.289. 

It may be noted that the pitching moment due to interference of the body upon the wing (—0.0848) is about equal 
in magnitude to the increment in pitching moment due to the normal component of the velocity induced by the 
wing upon the body (—0.0855), but the increment due to the longitudinal component (—0.1186) overshadows the 
other contributions. 


(E) CALCULATION OF THE OPTIMUM ANGLE OF WING Crw = (4/-/ M2 — 1)B2 + a6 


SETTING 
For the specific case under consideration we have 


(E-1) Derivation of Formula Giving Optimum Angle of 
Attack for the Wing ay 


—0.000207/0.1396 = —0.00148 


It may be assumed that a small deviation between +0.1659/0.1396 = +1.188 


the plane of the wing and the (x, y)-plane does not Pye coefficient of total lift, therefore, becomes: 
alter appreciably the equations given in the foregoing 


analysis. Thus it may be observed that it is possible to Cy = AiBi + Ape 


determine the optimum value for the angle of wing set- pore 

ting, i.e., the one that gives the minimum value of the 

total drag, for a given value of the lift coefficient, C,, of A; = (ACp/ Api) + a 

the wing-body system; however, the results should be A, = (4/V M,? —-1l) +a 


considered as merely indicative of the precise optimum. 


Let the drag coefficient of the body and that of the and the coefficient of total drag is: 


wing, both based on the same reference area, be denoted Co = Cop + An(B1:2) + Arz6i(B2) + Ace(B2?) 
by: 
h 
Con = Coo + Cow = BeCrw 
4 
respectively. The corresponding lift coefficients may 4,, = Ap =a+ta@, An = 
be expressed as: 4: M,?-1 
Crp = (ACz»/ABi) Bi + We may observe that 


and Bi = (C,/A1) — (A2/A1) Be 


A 
N- 
4 
roduced 
q 
‘ 
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and this expression, when inserted into the formula for 
‘Cp, produces: 


2 
Cp Cro + Au (a) | + 
| (a) | Bo + 
and thus 
dCp _ (8) 
A 
Aw | & (6 | — Ay (*) (B2) + Acs 


and when this is set equal to zero, there is obtained an 
expression for the optimum {2 which is: 


Gf 24 n(A2/A1) — As | 
A, 2A11(A2/A1)? 2A12 (A2/A1) + 2A 
Then, by examining the terms in d’?Cp/d{,’, it is 


found that d?C,/dB2* is > 0, because As. — Aj2 (A2/A1) 
is > 0, as may be seen from 


= 


2 — Ai A VM.? — 1 ({(ACp/ + ae] 
a,(ay + a2) 


(ACz,/ 481) + ae 


Consequently, the value of 62 given in Eq. (47) corre- 
sponds to an actual minimum of Cp. 


sag Optimum Angle of Wing Setting—Numerical Re- 
sults 


The optimum angle of wing setting is obtained, there- 
fore, as 


| 2Au(As/A1) — Av } 
2A1:(A2/A1)? — 2A12(A2/A1) + 2A» 
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If we assume that the slope of the lift curve* for the 
body is 2 (approximately), where the lift coefficient js 
based on the maximum cross-sectional area, then in the 
present case AC,,/ AG; = 0.561; and therefore 8, — 8, 
= 0.458C,. 

For instance, if C, = 0.3, then B; — Bo = 0.1375 
(= 7° 53’) and B, = 0.1522 (= 8° 43’), and, conse- 
quently, B. = 0.0147 (= 0° 50’). 

The optimum value of the angle of wing setting just 
found is probably somewhat high because of the sources 
of error already noted. The result just elicited, never. 
theless, must be regarded as a clear indication of how 
the body and the wing must be situated in order to re- 
duce the drag, for a given lift, on a specific wing-body 
system. 
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Problem Types in the Theory of 
Perfectly Plastic Materials” 


W. PRAGER? 


Brown Unwersity 


ABSTRACT 


Typical problems of the theory of perfectly plastic materials 
are discussed for the example given by the structure shown in 
Fig. 1. A simple graphical representation is developed which 
facilitates the study of the states of stress and residual stress pro- 
duced by slowly varying the load. 

For monotonically increasing load, three domains of mechanical 
behavior are defined: elastic deformation, contained plastic def- 
ormation, and unrestricted plastic flow. In the domain of con- 
tained plastic deformation the instantaneous stresses in the 
structure depend only on the instantaneous load and can be 
found from the minimum principle of Haar and von Karman. 

However, if unloading is permitted, the stresses depend on the 
complete history of loading, and the principle of Haar and von 
Karman must be teplaced by a minimum principle recently de- 
veloped by H. J. Greenberg. 

Cases in which the precise loading program is known before- 
hand are not often encountered in engineering. As a rule, only 
the extremes are known between which the load will vary. 
Under certain conditions, the structure will then ‘“‘shake down” 
toa state of residual stress such that all further variations of the 
load between the given extremes are supported in a purely elastic 
manner. This “‘state of residual stress’’ is independent of the 
precise program of loading. Conditions for the existence of such 
a state of residual stress are discussed, and a minimum principle 
is conjectured from which this state may be determined. 

Finally, structural stability in the plastic range is discussed. 
It is pointed out that the customary formulation of the stability 
problem for conservative systems (elastic range) is not adequate 
for nonconservative systems (plastic range). 


_ METHOD of theoretical or experimental stress 
analysis is based on a stress-strain law. Hooke’s 
law is generally accepted as an adequate basis for stress 
analysis in the elastic range. In the plastic range, 
however, no stress-strain law has found equally general 
acceptance. The selection of a stress-strain law ade- 
quate for the treatment of a specific problem is there- 
fore important. This selection is facilitated by a classi- 
fication of the problems encountered in the theory of 
plasticity. 

In the following, such a classification of problems is 
attempted for plastic materials that do not exhibit 
strain-hardening (perfectly plastic materials). To 
simplify the mathematical work as much as possible, 
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the various problem types will be discussed for the 
example given by the system shown in Fig. 1. The 
two-force members OA, OB, OC join the point O to the 
fixed points A, B, C. The system is symmetric with 
respect to the vertical OB, and the load P is acting 
along this axis of symmetry. The subscripts 1, 2, 3 
will be used to refer to the bars OA, OB, OC, respectively, 
and the forces acting in these bars will be denoted by 
Si, So, Ss. On account of the symmetry of the system, 
Si = S3. 

Consider, first, the mechanical behavior of the sys- 
tem in the elastic range. The elastic strain energy of 
the system is given by an expression of the form 


E = kS? + (1) 


where the positive constants k and K depend on the 
dimensions of the bars. A state of stress in the system 
(i.e., a set of values S,, S,) may be represented by the 
point with the rectangular coordinates 


x=SVk y= SVK (2) 


(see Fig. 2). This point will be called the “force point’’ 
and will be designated by F in Fig. 2. According to 
Eq. (1) the elastic strain energy of the system is then 
represented by the square of the distance of the force 
point from the origin. 

Equilibrium at O requires that S,x/2 + S, = P or 


xV/2/k + yV1/K = P (3) 


For a given value of P, this equation represents a 
straight line L in the x, y plane. This line will be called 
the ‘‘equilibrium line’; its distance from the origin de- 
pends on the intensity of the load P, but its slope is 
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independent of P. Any point on this line represents 
forces S; (= S3) and S:, which are in equilibrium with 
the given load. According to the principle of minimum 
strain energy, the state of stress actually produced in 
the elastic system by the given load corresponds to a 
smaller value of the strain energy than any other state 
satisfying the condition of equilibrium. Thus, the 
actual force point F is the foot of the perpendicular 
from the origin on the equilibrium line. If the load is 
increased, the equilibrium line moves away from the 
origin, maintaining its direction. Accordingly, the 
force point moves along a straight line through the 
origin. This line will be called the “line of elastic 
states.”’ 

Consider next the plastic behavior of the system. If 
the bars are made of a perfectly plastic material with 
equal yield limits in tension and compression, the abso- 
lute values of the forces S; (= S;) and S, cannot exceed 
certain yield values. If X and Y denote the corre- 
sponding values of x and y, the force point cannot move 
outside the ‘‘yield rectangle” with the sides x = +X, 
y = +Y. As the load is increased starting from zero, 
the force point moves along the line of elastic states 
until it reaches a side of the yield rectangle—e.g., the 
side x = X. In Fig. 2 the corresponding force point is 
denoted by F, and the equilibrium line by Z;. If the 
load is further increased, the bars 1 and 3 yield under 
constant stress. This means that the force point moves 
along the side x = X of the yield rectangle. When the 
equilibrium line has reached the position Z», for in- 
stance, the force point is at the intersection F2 of this 
equilibrium line with the line x = X. In spite of the 
fact that the bars OA and OC are now in the plastic 
state, their elongations remain restricted because the 
vertical displacement of the point O is determined by 
the elongation of the bar OB, which is still stressed 
elastically. Only when the force point reaches the 
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vertex F; of the yield rectangle will large plastic elonga- 
tions become possible. 

For monotonically increasing load, three domains of 
mechanical behavior must therefore be distinguished: 
elastic deformation (Fig. 2, segment QF, including end 
points); contained plastic deformation (segment F;F,, 
excluding end points); and unrestricted plastic flow 
(point F;). In the Theory of Structures only the first 
two domains are of interest, because large deformations 
cannot be tolerated. In the study of technological 
forming processes, on the other hand, only the third 
domain is of importance, because large deformations 
are desired here. 

Problems of contained plastic deformation have much 
in common with problems in elasticity, whereas prob- 
lems of unrestricted plastic flow are closer to problems 
in fluid dynamics. Accordingly, different theories of 
plasticity were developed in these two fields of applica- 
tion. The theories of plastic deformation (Haar and 
von Karman,' Hencky’) aim at problems of contained 
plastic deformation; the theories of plastic flow (Saint 
Venant,* von Mises,‘ Prandtl and Reuss*), at problems 
of unrestricted plastic flow. The stress-strain rela- 
tions of the former involve the strain itself; those of the 
latter, only the velocity strain. 

As Haar and von Karman! have pointed out, the 
principle of minimum strain energy of the Theory of 
Elasticity can easily be modified so as to apply to 
problems of contained plastic deformation. Indeed, the 
point /, in Fig. 2 is that point on the equilibrium line 
L, which is nearest to the origin without falling outside 
the yield rectangle. This agrees with the hypothesis 
of Haar and von Karman, according to which the actual 
state of stress corresponds to a smaller value of the elastic 
strain energy than any other equilibrium state of stress 
which satisfies the boundary conditions (prescribed sur- 

face stresses) and does not violate the yield condition. 


So far, it has been assumed that the load P increases 
monotonically. Suppose, now, that the load is de- 
creased after the state of stress represented by the 
point F, in Fig. 2 has been reached. For unloading of a 
two-force member, the decreases in load and elongation 
are related by Hooke’s law. Accordingly, unloading 
makes the force point move along a parallel to the line 
of elastic states. In particular, unloading from the 
state of stress represented by F: in Fig. 2 makes the 
force point move along FR. At the instant when the 
load is completely removed, the force point reaches the 
position R on the equilibrium line Ly through the origin. 
The point R therefore represents the residual forces 
after complete unloading from F2; it will be called the 
“residual force point.”’ 


From the preceding it follows that to each force 
point on the segment F,F; there corresponds a different 
residual force point. In other terms, the residual force 
point keeps moving along Lo as long as the permanent 
elongations of the bars S,; and S; vary but remains station- 
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ary whenever these permanent elongations remain un- 
changed. 

If unloading is admitted, the aforementioned modi- 
fied principle of minimum strain energy is no longer 
valid. It must be replaced by another minimum prin- 
ciple, the general form of which has recently been estab- 
lished by H. J. Greenberg (in a paper that is to appear 
in the Quarterly of Applied Mathematics). Applied to 
the example considered here and expressed in the ter- 
minology developed above, this principle takes the 
following form: As the equilibrium line undergoes a 
translation in accordance with the variation of the load, 
the force point describes the shortest path compatible with 
the conditions that it must remain on the equilibrium line 
and not leave the yield rectangle. Whereas the principle 
of Haar and von K4rméan establishes a direct relation 
between the instantaneous load and the corresponding 
stresses as they arise during the first loading of the 
structure, Greenberg’s principle establishes a relation 
between the rate of increase of the load and the corre- 
sponding rates of increase of the stresses in the structure. 
For plastic materials with strain-hardening, similar 
variational principles have been discussed by Prager® 
and, more recently, by Hodge and Prager.’ 

When unloading has occurred, the state of stress in 
the structure no longer depends on the instantaneous 
load alone but also depends on the complete history of 
loading. Unfortunately, cases in which the precise 
loading program is given beforehand are not often 
encountered in engineering. As a rule, only the ex- 
tremes between which the load will vary are known, 
but the precise manner of variation is not known. In 
Fig. 3, let L and L* be the lines of equilibrium corre- 
sponding to the extreme loads. If the structure is first 
loaded in tension up to the extreme load represented by 
the line Z, then unloaded and loaded in compression 
up to the extreme load represented by the line L*, the 
force point describes the path QFiF2F;F;. As the struc- 
ture is subjected to further cycles of loading, the force 
point traces the parallelogram F,F;F:F3. During each 
of these loading cycles plastic deformations occur while 
the force point describes the segments F;F, and F3F;. 

If, on the other hand, the extreme compressive load 
is represented by the line of equilibrium through F; 
(line L*’) the force point describes the path QF, FF; 
during the first loading and moves along F;F; from F; 
to F; and back to F; during further loading cycles. No 
plastic deformations occur during these cycles and the 
residual force point remains stationary at R. In other 
terms, the system ‘‘shakes down’’ to a state of residual 
stress such that all further variations of the load (be- 
tween the given extremes) are supported in a purely 
elastic manner. This state of residual stress to which 


the system will eventually shake down under loads 
that vary between given extremes is independent of 
the precise program of loading. 

General conditions for the occurrence of ‘shake 
down” in trusses have been discussed by Bleich* and 
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Melan,’ but the corresponding problem for other types 
of structures does not seem to have received much 
attention. For the example discussed above, Melan’s 
condition can be formulated as follows: The system will 
shake down tf there exists at least one parallel to the line of 
elastic states such that the segment which the extreme lines 
of equilibrium determine on this parallel lies within the 
yield rectangle. 

In an attempt at generalizing these results, one 
might, for instance, study the case of a cylindrical bar 
subjected to combined torsion and axial tension such 
that the torque and the axial force vary independently, 
each between given extreme values. The determina- 
tion of the shake-down state of residual stress in cases 
such as this constitutes an as yet unsolved problem of 
the mathematical theory of plasticity. If shake-down 
is at all possible for the given extreme loads, the state 
of residual stress to which the system shakes down prob- 
ably corresponds to a smaller value of the elastic strain 
energy than any other stite of residual stress which satis- 
fies the equilibrium conditions in the absence of loads and 
allows the body to support all further variations of the 
loads (between the given extremes) in a purely elastic 
manner. So far, this tentative theorem has not been 
proved, however. 

Finally, there is the problem of structural stability 
in the plastic range. The problem of column stability, 
for instance, is concerned with the determination of the 
smallest axial load for which the column can assume a bent, 
as well as a straight, equilibrium configuration. The 
classical treatment of column stability in the plastic 
range,'® however, is concerned with a more restricted 
problem—namely, the determination of the smallest 
axial load such that the column can pass from the straight 
to a bent equilibrium configuration while under the action 
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of this load. In the elastic range, the difference be- 
tween these two problems disappears, because the ex- 
terior forces necessary to maintain the system in a 
certain configuration depend only on this configuration 
and not on the preceding configurations of the system. 
In the plastic range, however, the history of deformation 
is of importance, and the two problems formulated 
above have different solutions."! 

These remarks can be illustrated on the example of 
the system shown in Fig. 1. Suppose that the con- 
figuration shown in Fig. 1 is assumed by the system 
under the influence of a compressive load P which 
already produces plastic deformations in the bars OA 
and OC but not yet in the bar OB. Then let this load 
P be increased to P + dP. If it is assumed that by this 
increase in load all three bars are further shortened, the 
forces in the bars OA and OC remain unchanged, since 
these bars yield in compression under constant loads. 
The new equilibrium configuration obtained under this 
assumption is symmetric with respect to the vertical 
OB. If, on the other hand, it {is assumed that the bar 
OA is lengthened while the bars OB and OC are short- 
ened by the increase in load, the force in OC remains 
unchanged, while the forces in OB and OC change by 
amounts proportional to the changes in length which 
these bars undergo. Further investigation reveals 
that, once the perfectly plastic bars OA and OC have 
reached the yield limit, the symmetric system of Fig. 1, 
subjected to a compressive load P, can absorb an incre- 
ment dP of the load in either of the two manners just 
discussed. If the symmetric configuration under the 
load P is denoted by J and the symmetric and unsym- 
metric configurations under the load P + dP by IJ 
and JI, respectively, the system can pass from J to 
either JJ or JI] provided that the load is increased from 
PtoP+dP. Asarule, the system can not pass from 


II to III under the constant load P + dP, however. 
For this to be possible, the load must considerably 
exceed the smallest load for which both a symmetric 
and an unsymmetric equilibrium configuration may 
first occur. 

Fig. 4 schematically represents this behavior of the 
system, as furnished by a linearized theory that neg- 
lects all but the first powers of the small displacements 
from the symmetric configuration J to the unsym- 
metric configuration //J. The ordinate in Fig. 4 repre- 
sents the load P; the abscissa, the lateral deflection yu 
of the point O. As the load P is increased, starting 
from zero, the equilibrium configuration remains sym- 
metric (i.e., the deflection « remains zero) as long as 
all three bars are stressed in the elastic range. When 
the bars OA and OC reach the yield limit (point a in 
Fig. 4), lateral deflection becomes possible under in- 
creasing load. The linearized theory leads to a linear 
relation between load and lateral deflection. In Fig. 4, 
this linear relation is represented by the sloping line 
through the point a. If the occurrence of lateral de- 
flection is prevented by guiding the point O along the 
vertical OB until the load has reached the value repre- 
sented by the point 6 in Fig. 4 and if the system is then 
allowed to assume an unsymmetric configuration, the 
relation between load and lateral deflection, as fur- 
nished by the linearized theory, is represented by the 
sloping line through b. The slope of this line is less 
than that of the line representing the bifurcation of 
equilibrium from a. Finally, for a certain load repre- 
sented by the point d, lateral deflection under constant 
load becomes possible. 

The linearized theory shows only that various equili- 
brium configurations are possible but does not indicate 
which of these will actually be assumed as the system 
is being loaded. It is to be expected, however, that a 
more complete theory would show that lateral deflec- 
tion will begin to develop at a and increase with in- 
creasing load following more or less the envelope of the 
sloping lines in Fig. 4. 

The kind of structural instability exhibited by this 
simple elastic-plastic system differs considerably from 
that encountered in completely elastic systems. This 
different behavior is caused by the nonconservative 
character of elastic-plastic systems. The point c in 
Fig. 4, for instance, can be reached from a along the 
sloping line ac. If lateral deflections are not permitted 
to occur at a and if the system is forced to remain 
symmetrical until the point 6 is reached, any lateral 
deflection will be accompanied by a further rise in 
load (sloping line through d), so that the point ¢ cannot 
be reached from 6. No general stability theory for 
nonconservative systems is available as yet, but the 
behavior of the simple system considered above is prob- 
ably typical for nonconservative systems in general. 

In conclusion, it seems worth while to point out that 
the preceding discussion applies in a qualitative manner 
also to the case of strain-hardening materials. Indeed 
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whereas the bars OA and OC of the simple structure 
discussed above were assumed to be perfectly plastic, 
the overall behavior of the structure, as represented by 
the relation between the load P and the vertical de- 
flection v of the point O, exhibits a sort of strain- 
hardening insofar as an increase of v requires an in- 
crease of P even when the bars OA and OC yield under 
constant stress. 
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Letter to 


Dear Sir: 

The following comments refer to the paper, ‘‘Range Perform- 
ance of Turbojet Airplanes,’’ by Irving L. Ashkenas (JOURNAL OF 
THE AERONAUTICAL SCIENCES, Vol. 15, No. 2, p. 97, February, 
1948). 

An important consideration that might well have been in- 
cluded in the paper is that of the relative values of the maxi- 
mum ranges which would be realized, under the assumed condi- 
tions, by cruising at the ceilings that correspond to the flight 
conditions defined by Eqs. (2), (6), and (10). Had such a com- 
parison been made, a welcome simplification of the practical 
prediction problem would have been indicated by the trivial 
differences between the ranges found for the three conditions. 
While this result may appear contrary to implications of the 
paper, it is entirely consistent with physical interpretation of the 
definitive equations. 

For convenience of reference, the three conditions of flight are 
identified as 


Case Condition Equation 
A Cp = 1.333 Copp (2) 

B Cp = 1.600 Cop (6) 

Cp = 1.500 Cop (10) 


In the paper of reference, it is shown that maximum range at 
fixed altitude is realized by maintaining the condition of 
Case A whereas, in the absence of the fixed height stipulation, 
maximum ranges are obtained in the troposphere and strato- 
sphere by cruising at the ceilings which correspond to Cases B 
and C, respectively. Let us consider the significance of these 
three flight conditions. 

Fig. 1 illustrates thrust-drag equilibrium for a jet-propelled 
airplane at the ceiling attainable in theeattitude for which D/ V— 
ie., T;/V—is minimum. This attitude is characterized by the 
telationship Cp = */;Cpp, the condition identified as Case A. 
Since constant values of dW/dR (lbs. of fuel consumed per 
mile flown) are represented by straight lines through the origin— 
defined by T = —(V/c)(dW/dR)—it is evident that maximum 
range is obtained at this, or any other, fixed altitude by flying 
at the lift coefficient which corresponds to Case A. 

Consider, now, the possibility of obtaining still greater range 
by cruising at heights above the ceiling for Case A. At the out- 
set, let us assume that the ceiling for Case A is above the tropo- 
pause; under this condition, the lines of dW/dR = const. re- 
main unchanged as altitude increases because ¢ is assumed to be 
constant in the stratosphere. Since an increase of altitude 


the Editor 


causes the 7, line to be lowered and the points of the 7; vs. V 
curve to be displaced horizontally, but nonuniformly, toward 
the right, it is evident that range will be augmented if, with in- 
creasing height, the intersection of the 7, line with the 7, curve 
moves into the region characterized by reduced values of 7/V 
and of dW/dR. 

Because the path traversed by the intersection point as alti- 
tude increases is not easily visualized in a Cartesian plot of T 
vs. V, the corresponding logarithmic chart, Fig. 2, is intro- 
duced to facilitate analysis. In this chart, all the lines of dW + 
dR = const. have the slope +1. Also, as height increases and 
density diminishes from o; to a2, the 7, line is lowered by the 
amount log (02/¢,), while the 7, vs. V curve is translated hori- 
zontally toward the right—without change of shape—by one- 
half that distance—i.e., the displacement is —'/2 log (o2/0,). 
Therefore, as Cx is increased—to maintain level flight at heights 
above the ceiling for Case A—the intersection point moves to 
lower values of 7/V (and of dW/dR) until it reaches the point 
of the 7; vs. V curve characterized by the slope +'/2. This 
slope corresponds exactly to the flight condition for Case C and 
the resulting improvement in range is 2.6 percent. Tocruise at still 
higher altitudes and lift coefficients would adversely affect range; 
for example, by going on to the condition for Case B, the gain 
would be reduced to 1.9 per cent. 

Treatment of the situation in which the ceiling for Case A is 
well below the tropopause involves an additional complica- 
tion: In Fig. 2, the lines of dW/dR = const. are displaced up- 
wardly as altitude increases, This results from the variation of 
c with o°-? in the tropopause. However, an analogous improve- 
ment of range is found but, in this instance, the maximum is ob- 
tained by cruising at the ceiling for Case B. In the troposphere, 
Case B yields 7.9 per cent more range than Case A. (In this re- 
gion, Case C yields only 7.1 per cent.) 


Fics. 1 and 2. 
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TABLE 1 
Relative Thrust Requirements, Ceiling Densities, and Ranges 
Basic Range Equation— 


dw c Cp ) 
A 


Case B 

Definitive equation (2) (6) (10) 

Co/Cop 1.333 1.600 1.500 

(0.333)'/2 (0.600) '/2 (0.500)'/2 
2.309 2.066 2.121 
(Ci! 0. 5699 0.5501 0.5606 

At Fixed Altitude— 

T,/T+A 1.0 0.8948 0.9186 

(dR/dW)/(dR/dW)a 1.0 0.9653 0.9837 

When Ceiling with W = W2 Is Below Tropopause— 

TAc/Te 1 1.172 1 

(dR/dW)./(dR/dW) ac 1 1.079 1.071 

When Ceiling with W = W, Is Above Tropopause— 

(dR/dW).'/(dR/dW) ac’ 1 1.019 1.026 


(Subscript c indicates ceiling.) 
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ANO RANGE RATIOS 


Fic. 3. 


The calculations which lead to the results cited above are 

summarized in Table 1. Fig. 3 presents, in another form, 

P graphical solutions of the typical problems that correspond 

a ies ; to heights above and below the tropopause. From this figure 
 &§ and the tabulated data, it is evident that the maximum ranges 
; attainable at the ceilings which correspond to the three flight condi- 
: tions differ by constant percentages so long as the tropopause is 
4 not crossed during cruising. [It may be of interest to note that 


, Fig. 3 reveals the cause of the ‘‘difficulty in determining the 
’ 35,000-ft. weight” from Eqs. (4) and (8)—to which the author 
4 refers on page 98. It is that immediately above the tropopause, 


maximum range is obtained by cruising at values of Cy, be- 
tween those for Cases B and C.] 


The foregoing comparison shows that, by deviating from the 
attitude in which 7/V is minimum (Case A), improvements in 
maximum range of 7.9 and 2.6 per cent might be attained in the 
troposphere and stratosphere, respectively—under the highly 
simplified conditions assumed. Lest the former figure be thought 
significant, consideration of the variation of dR/dW with height 
(see Fig. 3) will indicate the futility of attempting to obtain 
great turbojet range in the troposphere. At the higher levels, a 
calculated improvement of 2.6 per cent is certainly not significant 
when based on simplifying approximations of larger order. 

It thus appears that the rigorous analysis of the subject paper 
actually shows that the accuracy of turbojet range prediction is 
negligibly (if at all) improved by the consideration of flight condi- 
tions which differ from the one that yields maximum range at 
fixed height. It may be worth noting, specifically, that the 
equation 


23.8 T35 A W, 
R(mi.) = log 


C35 Cop 


which corresponds to cruising at the ceiling for Case A (above the 
tropopause) differs negligibly from the author’s Eq. (11) for 
Case C. Furthermore, even if such differences were not neg- 
ligible, recognition of the fixed ratios that, in stratosphere and 
troposphere, relate the maximum ranges for the three condi- 
tions of flight would have the welcome consequence of unifying 
the range criterion for variable height with that for fixed height. 
Therefore, as is evident from the basic range equation of Table 1, 
minimization of the ratio Cp/Cz,'/? is the fundamental aerody- 
namic criterion for the design of long-range turbojet aircraft. 

- In conclusion, the present writer wishes to commend the 
author for his clear demonstration of the compatibility of high 
speed with long range in turbojet aircraft, for pointing out (if 
not discovering) the constancy of optimum true cruising speed 
with diminishing weight in the stratosphere, and for calling at- 
tention to the remarkable improvement of range which occurs 
as the cruising altitude is increased. 

Before closing, attention should be drawn to a typographic 
error in the paper: On page 99, the value of (L/D),,,,. should 
read '/2VA,/Cpp. It is also noted that the numerical con- 
stant of Eq. (11) is based on the density at 35,000 ft. (0.000736) 
rather than on the value at the tropopause (0.000727) ; introduc- 
tion of this correction would make the constant 24.4. 


G. REID 
Professor of Aerodynamics 

Stanford University 

Stanford, California 
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Harmonic and Transient Motion of a 
Swept Wing in Supersonic Flow 


JOHN W. MILES* 
University of California, Los Angeles 


SUMMARY 


The results of previous papers on the harmonic! * and tran- 
sient? motion and loading of two-dimensional airfoils in supersonic 
flow are extended to include the infinite swept airfoil, provided 
that the component of flow normal to the leading edge remains 
supersonic. It is shown that the various normal force deriva- 
tives are reduced by the cosine of the sweepback angle and that 
the effective Mach Number corresponds to the component of 
flow normal to the leading edge. This result is related to the 
invariance of the magnitude of the circulation under sweepback 
and Busemann’s rule* with respect to compressibility effects. 
The flutter of a swept wing is considered, and it is pointed out 
that sweepback introduces additional aerodynamic coupling 
between torsion and bending. The flutter of a supersonic swept 
wing with subsonic leading edge is discussed, and an approximate 
treatment is suggested. 


NOTATIONS 
b = semichord measured parallel to free-stream velocity 
c = fraction of semichord of flap hinge line aft of mid- 
chord 
e = fraction of chord of elastic axis aft of line of aero- 
dynamic centers 
h = vertical deflection of wing measured at elastic axis, 


positive up 
j = (—1)'”, imaginary unit 
k = (wb/U), reduced frequency 
= distance aft of leading edge, measured parallel to 
free-stream velocity, cf. Eq. (10) 


t = time 

u = hyperbolic radius, cf. Eq. (2) 

w = vertical velocity in plane z = 0 

x, y,2 = coordinates (wind along x, y to starboard and z posi- 
tive up) 

2 = vertical displacement in ith degree of freedom, cf. 
Eqs. (21)—(23) 

Ai; |= dimensionless flutter coefficient giving force in ith 


degree of freedom due to motion in jth degree of 
freedom, cf. references 1 and 3 (see also end of this 
section) 

G(x,s) = Green’s function giving potential at x due to velocity 
at s, cf. Eq. (13) 

Jo = Bessel function, first kind, zero order 


L = lift, positive up 

M = Mach Number; also mid-chord stalling moment 

N = hinge moment, positive in sense to produce down flap 

S = area of disturbed flow intercepted by Mach cone 
from given point 

U = free-stream velocity 

a = angle of attack, measured as rotation about y axis 

B = (M? — 1)’; also flap deflection relative to airfoil, 
positive for down flap 

= (tano)/B 


Received November 14, 1947. 
* Assistant Professor of Engineering, Department of Engi- 
neering. 


—,n = source coordinates on (x, y) axes 
K = kM(M? — 1)7} 
= angle of sweepback of leading edge with respect to y 


axis 

= velocity potential 

0 = angle of twist about elastic axis (positive for leading 
edge up); also Mach coordinate, cf. Eq. (11) 

w = angular frequency 


The flutter derivatives Ars, in the absence of sweep, are defined 


by! 3 


where the a, refer to the following chordwise velocity distribu- 


s=h: 2,(x,t) = a,be 
S= a: ta(x,t) = —dabxet! 
s = B: t) = —agb(x — x>c 
= 
INTRODUCTION 


A ges NONSTEADY MOTION of an unswept two-dimen- 
sional airfoil has been considered in three previous 
papers.'~* The purpose of the present note is to extend 
the results of these papers to a swept wing whose lead- 
ing and trailing edges are ahead of the Mach waves— 
i.e., the case of pure supersonic flow. 


SOLUTION FOR A POINT SOURCE 


The solution to the wave equation for the linearized 
velocity potential of an oscillating point source moving 
through an ideal (nonviscous) fluid at a supersonic 
speed was found in reference 1, but the result, as stated 
therein, is in error. 

Following the analysis of reference 1, the possible 
solutions for a point source of strength g(t, )dt dn 
located at £, , 0 are given by 

do(x, y, 2) = —g(E, (1) 


u = [(x — §)? — — n)? — (2) 
Mk _ wb 


B= M*-1 (3) 
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where d is the unit of length and is taken as unity in the 
intermediate steps of the analysis. Eq. (1) is identical 
‘with Eq. (6a) of reference 1, except that both signs 
must be allowed u in the exponent, inasmuch as both 
solutions are permissible in the transformed coordi- 
nates in which the wave equation is solved. 

The vertical component of velocity due to the poten- 
tial of Eq. (1) vanishes everywhere in the plane z = 0, 
except at the source, where, following (analogously) 
Puckett’s analysis for steady flow,‘ it is given by 


w(x, 0,0) = — lim — g 
z—>0 O2/ a 


cos Ku + gp, Sin Ku 
E K K (4) 


Eq. (4) represents the integration of the source strength 
within the region cutout of the area dt dn by the Mach 
fore-cone subtended from the point (x,-0, 2). The 
source strengths of the two possible solutions given by 
Eq. (1) are specified by g, and g,, and the area inter- 
cepted by the Mach cone is bounded by the hyperbola 


B*(m? + 27) = (x — (5) 


(a is a constant, the exact value of which is immaterial.) 
Introducing the change of variable 


sin 0 = B(n’? + 27)'/2/(x — &’) (6) 


in Eq. (1) and carrying out the differentiation yields 
/2 


w(x, 0,0) = lim cos cos 0) + 
2+) /2 


Sin (xBz cos 0)|d0 = (7) 


Since the velocity potential must vanish with the 
velocity distribution w, it follows that the desired point 
source solution is 


y, 2) = w(t, cos xu (8) 
u 

[While Eqs. 16a and 18 of reference 1 give the exponen- 

tial in u, cos xu was implicitly used in Eq. (11), and all 


the subsequent results stated therein are correct.] 


SOLUTION FOR A LINE SOURCE 


The velocity potential at the point (x, 0, 0) on the 
swept, infinite span airfoil shown in Fig. 1* is given by 
integrating Eq. (1) over the Mach fore-cone subtended 
from (x, 0, 0): 


¢o(x, 0, 0) = we ne cos Ku ie 
S 
(9) 


where S is bounded by the Mach lines y = + [(x — £) + 
* In reference 1 the leading edge was placed at x = —1 instead 


of x = 0, but the latter choice is more convenient in the present 
analysis. 


Fic. 1. Plan view of wing swept ahead of Mach waves. 


8] and the leading edge = 7 tan ¢, if it is assumed that 
tan o < #—i.e., for the case where the component of 
flow normal to the leading edge remains supersonic, 
If it is further assumed that the flow is two-dimen- 
sional, and the transformations 


s=t-—ntane (10) 
cos = y + (1 — y’)B[n/(x — + ntano)] (11) 
y = tan (12) 


are introduced, Eq. (9) reduces to 
o(x) = —[1/8(1 — fo*w(s)G(x, s)ds (13) 


«M 
G(x, 5) = (as (x—s) cos x 
0 


cos | MEO 


The integral in Eq. (14) may be evaluated from the re- 
sult® 


54 cos (y sin = Jo[(x? + y*)!/2] (15) 
to yield 
. «M ran /2 
G(x, s) = e (x— s)(1 + B*y?)! | 
+) 
(16) 
From Eggs. (3) and (12), it is found that 
kM _ kM? _ _R(M cos o)? 
(l— (M?—1-—tan*oc) (Mcosc)?—-1 
(17a) 
kM sec o 
(1 — y?) (M? — 1 — tan? o) 
k(M cos (17b) 


[(M cos «)? — 1] 
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Introducing Eqs. (17) in Eqs. (13) and (16), it is seen 
that the dependence of the velocity potential on sweep- 
back is specified by 


o(x, k, M, = cos o¢(x, k, M cos a, 0) 


(17c) 


(18) 


APPLICATION OF BUSEMANN’S RULE 


It follows that the nonsteady motion of an infinite 
airfoil swept ahead of the Mach waves is equivalent to 
the motion of an airfoil having the same chord (2b) and 
chordwise velocity distribution [w(s)] im the direction of 
travel if the forces acting on the latter are reduced by a 
factor of cos o and its Mach Number is taken as M cos 
c, where M is the Mach Number of the actual flow. 
The various normal force derivatives, whether evalu- 
ated for a harmonic time dependence,' :* gust load- 
ing,? or an arbitrary time variation, are given by the 
same rule if the full velocity U is used in defining these 
coefficients, except where it enters the Mach Number. 

It should be particularly emphasized that the normal 
(to the leading edge) component of flow is used only in 
defining the effective Mach Number and not in defining 
k (since R = whb/U = wh cos a/U cos oa, the normal 
component of flow may be used in defining k if the chord 
is measured normal to the leading edge, and not in the 
direction of flow, as b is presently defined) or the dy- 
namic pressure. Hence, care must be exercised in 
applying Busemann’s® rule that the forces on a swept 
airfoil are determined only by the component of flow 
normal to the leading edge. If the chordwise slope of 
the airfoil is measured normal to the leading edge and 
no dihedral is present—i.e., if the angle of attack dis- 
tribution is made up of rotations about the swept axis— 
this rule may be interpreted literally, but if the angle of 
attack distribution is made up of rotations about an 
axis normal to the flow or if geometrical dihedral is 
present from any other source, the spanwise component 
of flow interacts with the dihedral angle to produce an 
additional upwash and contributes to the lift force act- 
ing on the airfoil. 

That the normal force derivatives with respect to an 
angle of attack measured as a rotation about the axis 
should be reduced by cos a, and not cos? a, is perbaps 
more easily seen by considering the normal (U cos a) 
and spanwise (U sin «) components of flow directly. 
Let it be assumed that the local pressure due to a normal 
component of flow V and upwash w is given by 


dF/dS = '/2pV?m(w/V) (19) 


Where m will, in general, be an arbitrary function of the 
geometry. Considering the swept wing, let a be the 
local angle of attack, measured as a rotation about the 
Now, if a is small, it may be regarded as being 


axis. 
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composed of two components, a cos ¢ and —a sin a, the 
former representing a rotation about the swept axis 
and the latter representing the effective dihedral with 
respect to the swept axis. The pressure due to the first 
component is given by Eq. (19), where V = U cosa 
and (w/V) = acoso, while the normal force due to the 
second component is given by Eq. (19), where V = U X 
cos gand w = (—asina)(—Usina); hence, the total 
pressure is given by 


iS = ; p(U cos a)*m E cos ¢ + 
(—asino)(—Usine)|_ 1 
| pU*(m cosa) a (20) 


Eq. (20) states that the effect of sweepback is to re- 
duce the local normal force derivative by cos ¢. Its 
derivation does not replace the derivation of Eq. (18), 
inasmuch as it is assumed that m itself is invariant with 
respect to sweepback, whereas, in actuality, it too must 
be modified in accordance with Busemann’s rule when 
compressibility is considered. 

Eq. (20) may also be interpreted as stating that the 
magnitude of the circulation about the swept axis is 
unchanged by sweepback, the cos o entering from the 
vector cross product in the Kutta-Joukowsky law 
(dF/dS = pU X ¥ where y¥ is the local circulation). 
While such a postulate (that only the normal component 
of flow need be considered insofar as compressibility 
considerations are concerned) might be used to state 
Eq. (18) directly, particularly in retrospect, this is per- 
haps asking a little too much of physical intuition. 
Nevertheless, both rules appear to be generally appli- 
cable to problems of both steady and nonsteady flow in 
the case of purely subsonic or purely supersonic lin- 
earized two-dimensional flow and to many cases of 
three-dimensional flow; in the case of mixed supersonic 
flow, the rules may be locally applicable, but the varia- 
tion of m may be complex. 


FLUTTER CONSIDERATIONS 


In treating the flutter of a swept wing, it will be 
generally convenient to describe the structural defor- 
mation in terms of bending (#) about a section perpen- 
dicular to the swept axis, twisting (@) about the swept 
elastic axis, and control surface rotation (8) about a 
hinge line parallel to the swept axis; the deflections 
with respect to an axis perpendicular to the free-stream 
velocity may then be written (neglecting taper effects 
and assuming all angles to be small, so that they may be 
treated as vectors) 


y’,t) = | + 2e cos ay) | be**' (21) 
z.(s,y’,t) = — cos — sine | x 
oy 
b(s —1)e* (22) 


(12) 
(13) 
(14) 
he re- 
(15) 
(16) 


za(s, y’, t) = — cos oB(y’)b(s — 1 — 
s>(1+o) (23) 
s<(1+6) 
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where y’ is measured along the swept axis, s is measured 
in the direction of flow [cf. Eq. (10)], and e, the eccen- 
tricity, is the fraction of the chord which the elastic 
axis lies behind the line of aerodynamic centers. 


It then follows the definition of the aerodynamic flutter derivatives, in Eqs. (22) and (23) of reference 1 and 
Eqs. (4)—(6) of reference 3,* that the forces and moments on a swept supersonic airfoil, based on strip theory, are 


val + AcB(k, M cos cos (24) 


given by 
dL(t) cos =| { [* ] 
= Ach(k, M —+2 Ay’ 
dy ch( cos + 2e cos oA(y’) | + 
Aca(k, M cos cos — 
dM(t) _ 


cos 


dy B(M cos oc) 


dant M cos (y') + 2e cos 


Aaa(k, M cos a] cos — sing + AaB(k, M cos a) cos (25) 


dy B(M cos oc) 


Aba(k, M cos cos — sing 


dN(t) cos Abhi, M cos (y’) + 2¢e cos | + 


) 


+ AbB(k, M cos c) cos (26) 


where (dL /dy) is the lift per unit span, (d//dy) is the mid-chord stalling moment per unit span, and (dN/dy) is 


the hingé moment per unit span, y being measured perpendicular to the free air stream. 


B = (M? cos? « —1)'”. 


B (M cos a) indicates 


The results of Eqs. (24)-(26) are also applicable to a purely subsonic swept wing, provided that the appropriate 
flutter derivatives” ' are used and 8 (which will be imaginary for M cos ¢ < 1) is replaced by its absolute value. 
It is evident from Eqs. (24)—(26) that sweepback introduces an aerodynamic coupling between torsion and bend- 
ing even when the elastic axis coincides with the line of aerodynamic centers. 


* See section on notation for definitions in absence of sweepback. 


SuUBSONIC LEADING AND TRAILING EDGES 


The general result of Eq. (9) is valid for any airfoil 
in the supersonic régime where a linearized treatment 
of the flow is allowable; however, in those cases where 
either the leading or trailing edges of the planform are 
swept inside the Mach cone—i.e., where the component 
of velocity normal to these boundaries is subsonic 
(generally called mixed supersonic flow)—the area S 
includes disturbed regions off the wing (analogous to 
the wake region in purely subsonic flow) over which the 
velocity distribution is not explicitly known, and, in the 
case of the subsonic trailing edge, _ Kutta condition 
must be satisfied. 


While the problem of the two-dimensional oscillating 
airfoil in mixed supersonic flow appears to be difficult 
and has not, to the author’s knowledge, yet been solved, 
the solution of particular steady flow problems’ indi- 
cates that the flow should be approximately subsonic 
over most of the span of an actual wing. Accordingly, 
the results of Eqs. (24)-(26) could be used for a wing 


swept behind the Mach cone, with the supersonic 
derivatives': *.* being used over the inboard position 
of the span and the subsonic derivatives (at J = 0, 
or preferably at ./ cos o) over the outer portion of 
the span. The dividing section between the two types 
of flow is probably best determined by an examination 
of the pressure distribution over the same wing in 
steady flow. In following such a procedure, it should 
be recalled that those supersonic flutter derivatives in 
the pitching degree of freedom are referred to the mid- 
chord, whereas the corresponding subsonic derivatives 
are referred to the quarter chord. 
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BuCKling of a Column 
elastic 

T. H. LIN* 
ce 1 and 

COTY, are 


ABSTRACT 


In this paper, the buckling of a column with equally spaced 
elastic supports under compressive load is investigated, and an 
equation is derived for any number of equally spaced elastic 

supports. It is assumed that the elastic supports give only 
} (24) linear but not rotational restraints on the bar and the bar itself 
is in long column failure. As the spring constant of the elastic 
supports increases, waves of buckling increase until there is one 
half wave between every two adjacent supports and until further 
increase in the spring constant does not increase the buckling 

load. This constant is called the maximum critical elastic 
\ (25) constant. An equation for calculating this constant is given. 

Curves of the ratio of buckling load to Euler load vs. spring 
constant for different modes of buckling are shown so that 
with a given column on equally spaced elastic supports of known 
spring constant, the buckling load can be readily calculated. 


(26) NOTATIONS 
l = span between supports 
V/dy) 1S EI = flexural rigidity of the column 
idicates P = end load 
Pr = Euler’s load per span = 7r°EJ//? 
. P = critical buckling load 
ae) 
pra K = elastic constant of elastic supports 
value. K, = critical elastic constant of elastic supports 
d bend- Wi = total span of the column with elastic supports 
Ve = energy due to bending 
Ve = work done by elastic supports 
Vp = work done by load P 
p = n/W, a parameter 
m,n,q = constants positive integers 
ersonic THE GENERAL SOLUTION OF 
osition MULTIELASTIC SUPPORTS 
= 0; 
ion A COLUMN BAR with multielastic supports 
> types in — axial load “6 shown in Fig. 1. Assuming 
ination e buckling curve to be 
ying in 
= ZA, sin (mrx/W1) (1) 
e mid- | Which satisfies the boundary conditions at x = 0 and 
rations x = WI. Then 


4 wil 2 
Ve = f ( A,n? sin dx (2’) 
2W'l* Jo Wl 
illating 2 Wi - -\2 
ciences, Vp = ( A,Mn Cos dx (3’) 
0 . Wl 


re Using the principle of orthogonality, after integrating, 
ilating Received November 14, 1947. 


*S.M., A.F.R.Ae.S. Lieutenant Colonel. 
7 BS. Captain. 
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with Elastic Supports 


K. S. CHING? 


Chinese Air Force 


= (2) 
n=1 
= (3) 
n=1 


Let m be a positive integer from 1 to IV; then, 


Ww 
if K = = mil 


26 +4) (4) 


m=1 n=1 


or 


By Rayleigh’s principle, 

Vp = Vk + Vp (5) 
Substituting the values of Vz, Vp, and I’x from Eqs. 
(2), (3), and (4) into Eq. (5), we get 


2 4 


1K Ms sin (6) 
m=1 n=1 


Introducing S = 7°*/2WIK rt solving Eq. (6) for P., 


n=1 
According to the theory of minimum energy, the 
constant A, should have such values that make P, a 
minimum value—that is, 


P, = 


OP./0A, = 0 
where g = 1, 2,3, ..., ©. The result is 
Pe) 
Sq’A sin 


n=1 


Fic. 1 


Vp 
V 
P 
= 
n=1 n=1 
4 
x 
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Define 


= — (8’) 


K = C(P;/l) 
where C is a constant. Then, 


S = (x?/2WC)(1/Pz) 


Therefore, 


2 2 P 
owe \w P, (8) 


As the value of g varies from one to infinity, there js 
an infinite number of linear equations for determining 
the value of A,—i.e., 


bide + = 0 


where g = i, 2, 3, ..., &. 


In order that the coefficients A, cannot equal to zero simultaneously, the determinant of the group of linear 
equations must be zero. Furthermore, since the sine series is used to express the deflection curve, (2W — 1) 
terms of the sine series would give a result accurate enough for practical use. The error caused by limiting the 
infinite determinant arbitrarily to the first (2W — 1) terms will be estimated by comparing the results so obtained 
with exact results as shown in reference 1. The determinant is: 


¥(2, 1) + ¥(2, 2) 
¥(2W — 2,1) ¥(2W — 2, 2) 
¥(1,2W — 3) v¥(1, 2W — 2) ¥(1,2W — 1) 
¥(2,2W — 3) y¥(2, 2W — 2) ¥(2,2W — 1) 
¥(3, 2W — 2) ¥(3, 2W — 1) 
WOW = 3,2W-2) (10) 
¥(2W — 2,2W — 3) saa 3 = ¥(2W — 2,2W — 1) 
v2W-1,2W-3) ¥2W-1,2W-2) 
2W — 1) 
It is easily seen in the determinant that all of the 1 m iii 
y functions can be classified into four groups: (A, B) = : ) cos (A — B) a 
=¥(A = B= W) iv 
= ) cos (A + B) — 
A+B... 2 W 
+ integer 
3,48 ) since 
Furthermore, it can easily be seen from the determinant ) cos n0 = — 1, 1sin (2n + 1)(6/2) 
that all the y functions associated with ¢ functions in Some 2 2 sin 6/2 
the principal diagonal are in the form of yj, and all 
y functions in the secondary diagonals are in the form therefore, 
of ¥4, with y¥2 in the center of the determinant. 
All of the other elements in the determinant are in W(A, B) = */s(a — 7) 
the form y3. 
By the definition of the function of y, where 
w sin (2W + 1)[(A — B)x/2W] 
B) = > “sin gin, sin [(A — B)x/2W] 
fom W W 


Therefore, by trigonometry, 


_ sin (2W + 1)[(A + B)x/2W] 
sin [(A + B)x/2W] 


Fo 


Fot 


The 

of tl 
That 
elast 


gives 
elast 


Th 


and 


Ama . Bmr 
=e ae W(A, B) = ) sin —— sin —— (7) 
m=1 
and 
Fo 
Let 
| 
tie, 
; 
At 
ing 
whe 
struc 
begi: 
(W/ 
R 
4 
Ai, 2 
book 


(8) 


y, there js 
termining 


of linear 
— 1) 
iting the 
obtained 


0 (10) 


W 


Im 


BUCKLING OF A COLUMN WITH ELASTIC SUPPORTS 


For Vy 
a+2W+1, 7r=1, vi = W/2 
For ¥2, 
a= 2W+1, 7=2W+1, = 0 
For 
¢: + (W/2) 0 0 0 
0 ¢: + (W/2) 0 0 
0 0 aha 
ow 
0 0 
0 —-W/2 0 0 
—W/2 0 0 0 


In order to expand this determinant, we define 
Ai, =w—1 aS the above determinant; then, by eliminat- 
ing outer rows and columns, we get 


Ayew-1 = + + w)_(#) | ow —2 


where Ao, 2xv—2 is the determinant with the same con- 
struction as in Aj, 2-1 but where the principal diagonal 
begins with ¢. + (W/2) and ends in the form of ¢2y—2 + 
(W/2). 

Repeating this process, we get, finally, 


The function ¢y = 0 gives the buckling load with each 
of the span of the whole bar as a Euler's column. 
That means there are no lateral movements of any 


elastic supports. The part 

Tl W W\ (W\]_ 

on dew —n + 9 7) = 0 
gives the buckling loads for various deflections of the 
elastic supports. 
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a=(-1)At4, += (-1)At = 0 
For y¥4, 
a=1, r= 2W +1, = —W/2 


Therefore the determinant may be simplified by intro- 
ducing the values of y as follows: 


0 —W/2 
—W/2 0 
0 0 
pene = 0 (11) 
0 0 
daw -2 + (W/2) 0 
0 dow -1 + (W/2) 


The above equation may be written in the form 


Ww-1 
+ (W/2) (bn + = 0 


n=1 


(12) 


where the function ¢ is defined by Eq. (8). 

Eq. (12) is the final result of this buckling problem; 
the number of complete buckling waves is n/2, and, 
consequently, the wave length is 2W1//n. 


VALUE OF MaximuM CRITICAL ELASTIC CONSTANTS 


If the value of C,,,,. is the maximum critical value 
of the elastic constant the bar will buckle in the form 
that has maximum number of waves. Therefore, 


Cnaz. may be obtained by putting P, = P, and n = 
W — 1in Eq. (12). We get: 


— 9W! + 6W? 
5we + Ws 


and the value of C,,,,, when W approaches infinity is 


(13) 


lim = = 3.95 


The exact value is ‘‘4,”’ as discussed in reference (1); 
the error is nearly 1 per cent. 


COMPARISON OF 8 AND C 


The 8 in Timoshenko’s book equals 1/C in the present paper, so the following table is a comparison. 


W 2 3 4 5) 6 7 9 11 
8 0900 0333 0293 0276 0268 0263 0258 0255 
1 

C 0904 0337 0296 279 0271 0266 0261 0259 


and the result C = 2X[1 + cos (x/IW)] by Ebner’s formula gives the value of 1/8 as tabulated in Timoshenko’s 
book.! 


B) mr ; 
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The values for C,,,,, up to 20 elastic supports have been calculated with Eq. (13) and plotted in Fig. 2. 


NUMERICAL CALCULATIONS 


Take an element out of Eq. (12) and let p = 


n/W. Then, after manipulation, 


[o* — (P./Pz)[(2 — p)? — (P./Px)] 


C = — p)z]? 


+ (2 — — + (2 — 


C wax. / 


Fic. 2. 


The difficulty is to find what is the value of p for 
obtaining the minimum value of P,/Pz. Evidently 
an envelope is desirable for this purpose, but it is 
difficult to obtain such an equation. Therefore, a 
set of curves for p varying from 0.1 up to 0.9, which 
covers all the range of p since the range of m is from 
1 to W — 1, has been computed and is shown in Fig. 
3. This figure will help the stress analyst to find the 
value of p for a particular value of W, and with a few 
trials it is easy to get the correct answer by using 
Eq. (14). 


Example 
Given Data: 
I = 0.004 in.‘ 
1 = 21 in. 
E = 10 X 106 lbs. per sq. in. 
K = 75 Ib. per in. 
W =6 
A = 0.04 sq.in. 
10’ 0.004 


21 X 21 


os pire 
Pc 
Pe Va 
1 2 3 
Fic. 3 
900 75 
— = — = 427 lbs. C 
i 7 Ibs. per in., 27 1.75 


= 0.83, p = 


The values of p in this case are: p = 5/6 
= 0.33, and p = 


4/6 = 0.67, p = 3/6 = 0.5, p = 2/6 
1/6 = 0. 17. 

From Fig. (4) we see that the minimum value of p 
is about 0.6. Therefore, use p = 0.67. Substituting 
into Eq. (12) or Eq. (14) or estimating from Fig. 4. 


P./Pz = 0.80 


Therefore, P, = 0.80 K 900 = 720 lbs. 
The corresponding stress 


S = 720/0.04 = 18,000 Ibs. per sq.in. 
The number of complete buckling waves is 4/2 = 2. 
The wave length = (2)(6)(21)/4 = 68 in. 
REFERENCES 


! Timoshenko, Theory of Elastic Stability. 
2 Aitkin, A. C., Determinants and Matrices. 
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Machmeters for High-Speed Flight Research 


HERBERT SCHAEFER* 
Kollsman Instrument Division of Square D Company 


SUMMARY 


The Machmeter, an instrument designed to indicate the Mach 
Number at which an aircraft is flying, is being used as a test 
instrument as well as a flight instrument. The objective of this 
paper is to give a résumé of the theory and design of various types 
of Machmeters now in use. 

After an explanation of the Mach Number and a brief survey 
of transonic and supersonic flow conditions, the aerodynamics of 
the Machmeter is presented. For the supersonic speed, a sim- 
plified exact equation is given for the stagnation pressure as a 
function of the flight Mach Number and the free-stream pressure. 
Graphs are presented showing the relations of pressure and tem- 
perature at the stagnation point vs. the flight Mach Number. 

After a short description of the design of a subsonic Mach- 
meter, a more detailed description of the newly developed tran- 
sonic Machmeter follows, mentioning a number of design details. 
Then the Mach Number Limit Switch is briefly discussed, fol- 
lowed by a description of the Maximum Allowable Airspeed or 
Mach Airspeed Indicator, the aerodynamic theory of which is 
also presented. 

The location of the pressure vents is discussed, and a new type 
of supersonic static pressure probe is suggested. 

The possibility of a Machmeter operated by temperature pa- 
rameters is mentioned, but it is considered not feasible because 
of the difficulties in obtaining the proper temperature values in 
flight. 

With the experience already gained in developing Machmeters, 
it seems certain that, in order to keep step with the pace of the 
increase in flying speeds, Machmeters with increased ranges of 
Mach Number and altitude can be designed and built. 


(1) INTRODUCTION 


_ INCREASE in flying speed during the last decade 
has brought into practical aerodynamics as a new 
parameter, the Mach Number, the ratio of true air speed 
to the speed of sound. The dependence of aerodynamic 
phenomena on the Mach Number has been known to 
aeronautical engineers for a long time, but for the flying 
speeds attained until recently, this dependence on the 
Mach Number was negligible, and so, except in the 
domain of ballistics, the Mach Number was chiefly of 
theoretical interest. However, as flying speeds reached 
such magnitudes that in certain regions around an air- 
craft local air speeds were attained which approached 
or even exceeded the local speed of sound, it became im- 
perative to know the Mach Number at which the air- 
craft was flying in order to determine the safe opera- 
tional speed limits. For this purpose the Machmeter, 
an instrument designed to indicate the Mach Number 
at which an aircraft is flying, was developed. For 
operational use in high-speed aircraft, a sister instru- 


Presented at the Instruments Session, Sixteenth Annual 
Meeting, I.A.S., New York, January 26-29, 1948. 
* Project Engineer. 


ment of the Machmeter, the Maximum Allowable 
Airspeed Indicator or Mach Airspeed Indicator, was 
also developed. The theory and design of various types 
of Machmeters is the topic of this paper. 


(II) CoNcEPT OF THE Mach NUMBER 


The Mach Number is defined as the ratio of the flow 
velocity of a fluid to the velocity of sound in that 
fluid. It was named in honor of the Austrian physicist 
Ernst Mach, who worked on basic research in high- 
speed aerodynamics and ballistics as early as the last 
century. In the case of an aircraft, it is the so-called 
“Free-Stream Mach Number” or “Flight Mach Num- 
ber’ that is of chief interest; its definition is, the ratio of 
the true air-speed V of the aircraft with respect to 
the undisturbed air ahead to the speed of sound in 
the same air c: 


(1) 


The term true air speed is self-explanatory. The speed 
of sound in air is equal to the velocity of propagation of 
small pressure disturbances through the air. 

Thermodynamic changes occurring during the propa- 
gation of such small pressure disturbances can be con- 
sidered as being adiabatic, since all the changes of the 
parameters involved are small and no heat exchange 
with the surrounding air occurs. Assuming, further- 
more, that the air is a perfect gas, which follows the 
general thermodynamic equation of state p/p = RT, 
the speed of sound is given by the equation 


¢ = Vdp/dp = Vkp/p = VkRT (2) 


In these equations the symbols represent the following 
thermodynamic parameters: p, the barometric pres- 
sure; p, the mass density; 7, the absolute tempera- 
ture; R, the gas constant; and k = c,/c,, the ratio of 
the specific heats at constant pressure and constant 
volume, the numerical value of which is usually taken 
for air to be k = 1.400. 

Since in the atmosphere, pressure, density, and tem- 
perature are all variables, the speed of sound is a variable 
also. As the temperature diminishes with altitude, the 
speed of sound also diminishes, since it varies as the 
square of the absolute temperature, as indicated by 
Eq. (2). Under standard conditions, the magnitude 
of the velocity of sound is 761.06 m.p.h. at sea level 
and is 662.23 m.p.h. at the tropopause—i.e., an alti- 
tude of 35,332 ft. This value remains approximately 
constant up to 100,000 ft. In the upper atmosphere 
the temperature and, consequently, the speed of sound 


M = V/c 
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will increase and decrease several times with increase 
in altitude, as recent research has shown.! 

From the dependence of the sonic speed on altitude, 
it is evident that two aircraft flying at a certain true 
air speed at two different altitudes will fly at different 
Mach Numbers. If an aircraft is diving with constant 
true air speed, the flight Mach Number will continu- 
ously decrease because the speed of sound is increasing. 
This fact has probably saved the life of many a pilot 
who was diving, with his aircraft out of control, because 
it was flying faster than its critical Mach Number; 
during the dive the flight Mach Number became sub- 
critical, and so the pilot could regain control of his 
aircraft. 

The physical significance of the Mach Number be- 
comes evident if Eq. (1) is squared: 


M? = V2/c? = V2/kRT 


for then, the numerator represents a direct measure of 
the controlled or directed kinetic energy of the flow. 
The denominator represents the random translatory 
and rotary kinetic energy of the molecules, because it 
contains the absolute temperature which measures that 
energy. The square of the Mach Number represents, 
therefore, the ratio of the directed kinetic energy to the 
random molecular energy.* Another manner of de- 
fining the square of the Mach Number is that it repre- 
sents the ratio of the mass forces to the elastic forces.® 
It is interesting to note that in aerodynamic equations 
the Mach Number rarely occurs in the first power but 
is almost always squared. 


(III) FLow PHENOMENA IN TRANSONIC AND 
SUPERSONIC FLow 


A little digression into the fundamentals of high- 
speed aerodynamics will show at once the great impor- 
tance of the Machmeter as a tool for high-speed flight 
research. This is particularly true for transonic speeds, 
in which the flow around the aircraft is a mixture of 
subsonic and supersonic flow. Since both these types 
of flow follow different physical laws, a particularly 
difficult flow pattern results. It is chiefly for the re- 
search in this transonic region that a transonic Mach- 
meter with a range of flight Mach Numbers 
0.5 to 1.5 was developed. 

If an airflow takes place at extremely slow speed— 
that is, at a low Mach Number with subsonic flow 
velocity—the air can be assumed to be incompressible. 
Up toa Mach Number of 0.15—.e., an air speed of only 
15 per cent of the speed of sound—the errors introduced 
by the assumption of incompressibility are of a mag- 
nitude of about 1 per cent. Since an incompressible 
fluid has no elasticity, the influence of the Mach Num- 
ber, which represents the ratio of the mass forces to the 
elastic forces, is negligible. The important parameter 
for such flows is the Reynolds Number, the ratio of 
the inertia forces to the viscous forces. If the speed 


of flow is increased and, therefore, the Mach Number 
is also increased, the effect of compressibility becomes 
increasingly greater. When the flow velocity is jn. 
creased still more, the Mach Number finally becomes 
greater than unity—that is, the flow velocity is now 
greater than the velocity of sound and is therefore 
supersonic. Flow phenomena that occur in this domain 
are entirely different from those in subsonic flow. The 
most important of these are the following: 

Pressure disturbances cannot be propagated into 
every part of the fluid but are restricted to certain do- 
mains within the fluid which lie downstream from the 
point of disturbance. 

Flow through convergent stream tubes is decelerated, 
while flow through divergent stream tubes is acceler. 
ated; this is exactly the opposite of the behavior of 
subsonic flows. 

Furthermore, discontinuities with an increase in 
entropy, or so-called shocks, can occur in supersonic 
flow; these are impossible in subsonic flow. 

The reasons for this entirely different behavior of 
supersonic flow is, mathematically expressed, that 
the differential equations of motion for supersonic flow 
are hyperbolic, while those for subsonic flow are 
elliptic.‘ 

It is particularly the phenomenon of the occurrence of 
shocks that causes the most difficulties in high-speed 
flight. A shock is actually nothing but a special de- 
tached interior boundary layer in a fluid, which repre- 
sents a discontinuity sheet through which the pressure, 
density, temperature, and entropy increase suddenly 
while the flow velocity drops.’ Such shock dis- 
continuity sheets are extremely thin. Taylor and Mac- 
coll, for instance, give the thickness of a shock wave 
through which the velocity drop is approximately '/y» 
of the speed of sound as less than 0.0033 mm. = 0.00013 
in.© Such shock waves will occur as soon as the local 
sonic speed is reached somewhere near the aircraft. 
Usually they will appear first on the upper surface of 
the wing and then, as the speed is increased, on the 
lower surface as well. The reason for their occurrence 
lies in the fact that pressure disturbances from the 
trailing edge of the airfoil cannot pass through the 
region of supersonic flow, because the speed of their 
propagation is only that of the speed of sound. The 
pressure disturbances, therefore, accumulate at the 
boundary between the subsonic and the supersonic 
flow and form a shock wave. Such shock waves have 
the tendency to distort the entire flow pattern around 
an aircraft and simultaneously disturb the pressure dis- 
tribution around the wings. The creation of a region of 
high pressure behind the shock will cause adverse pres- 
sure gradients within the boundary layer, because 
a propagation of pressure through the subsonic part 
of the boundary layer from the downstream to the up- 
stream side of the shock is possible. This phenomenon 
usually causes boundary separation. The occurrence 
of shock waves on the lifting surfaces causes, in gen- 
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eral: decrease in lift; increase in drag; pronounced 
changes in aerodynamic stability; loss of effectiveness 
of the control surfaces, which in extreme cases may 
even reverse their intended effects; structural dam- 
ages to the control surfaces and other parts of the air- 
craft because of violent oscillations due to shock stall; 
and severe buffeting.””* 

The Mach Number of the undisturbed flow, or 
fight Mach Number, at which the maximum local 
Mach Number becomes unity in some region around 
the aircraft, is called the lower critical Mach Number. 
The disturbing phenomena listed above will not 
occur immediately at this particular Mach Number, 
but at one somewhat higher than the critical flight 
Mach Number, although still below unity.’ How- 
ever, this critical Mach Number is not valid for all 
fight conditions; the aircraft will have several critical 
Mach Numbers: one for straight and level flight, one 
for diving, and soon. The lower critical Mach Number 
isan important parameter for the pilot of a high-speed 
aircraft, because it warns him that he entered the tran- 
sonic range with its impending dangers. 

At a certain flight Mach Number the flow around an 
aircraft becomes entirely supersonic, except for the 
boundary layer and small regions near the nose of the 
aircraft and the leading edges of the airfoils. This 
fight Mach Number is called the upper critical Mach 
Number. This parameter is also important to the 
pilot, because it indicates the end of the transonic range 
and marks the beginning of completely supersonic flow 
around the aircraft.*® 

The lower critical Mach Number of an aircraft and 
many of the detrimental compressibility phenomena 
can be ascertained with the subsonic Machmeter, 
which has a maximum range up to Mach Number 1.0. 
This instrument has already been in use for many 
years. It is now up to the newly developed transonic 
Machmeter, the range of which is extended beyond 
Mach Number 1.0, to assist the test pilot in forcing 
his way beyond the transonic flight region into the 
domain of completely supersonic flow, which until 
today has been the realm of meteors, jet-propelled 
missiles, and rockets. 

Machmeters are not exclusively flight research in- 
struments; they are also used operationally. Further- 
more, they have proved their usefulness for aerody- 
namic research in wind tunnels. 


(IV) AERODYNAMICS OF THE MACHMETER 


There are several possibilities for obtaining the free- 
stream or flight Mach Number by measuring various 
aerodynamic parameters in flight. One method makes 
use of the stagnation pressure 3 at the head of a pitot 
tube and the static atmospheric pressure /;. These 
two parameters may be used either as differential 
pressure Ap = p; — pi, as in the case of the conventional 
air-speed indicator, or as a pressure ratio p3/p;. An- 
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other possibility is to use the stagnation temperature 
T; and the free-stream temperature 7}. 

In the following section the aerodynamic relations 
between these parameters and the flight Mach Number 
will be developed. Since the subsonic and the super- 
sonic flow conditions are different, both cases are 
treated separately. 


(A) Subsonic Air Speed 


In the mathematical treatment of this problem the 
following simplifying assumptions can be made: 

(a) The air is considered to be a continuous medium. 

(b) The air is considered to be an ideal gas obeying 

the general thermodynamic equation of state. 

(c) Being an ideal gas, as assumed under (b), its 
viscosity is zero. 
The airflow is considered one-dimensional; 
this assumption is justifiable, because only a 
small stream tube is being considered, which is 
leading directly to the small orifice at the head 
of the pitot tube. 
The coordinate system of reference is fixed in the 
aircraft, which is traveling with constant velocity 
through the undisturbed air; this defines the 
resulting flow as a steady flow, and all time 
derivatives vanish. 
(f) The motion is considered irrotational. 
(g) No outside forces, such as gravity, are influenc- 

ing the flow. 

(h) The flow occurs without any exchange of heat 
with the surroundings—i.e., it is adiabatic. 

(i) The entropy of the air remains constant—i.e., 

the flow is isentropic. 

With the above assumptions the general equation of 
energy in the following form can be used as the starting 
point of the derivation.” 


(d) 


(e) 


Schematic diagram of a pitot tube at A, subsonic, and B 
supersonic speed. 
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Fic. 2. Bow wave ahead of a pitot tube within a \-shock in a 
wind tunnel at M = 1.225. (Fig. 4 of reference 10.) 


u(du) + di = 0 (3) 


In this equation u represents the velocity in the flight 
direction of the aircraft and 7 the heat content or en- 
thalpy of the air. By designating the free-stream con- 
ditions with the suffix 1 and the conditions at the 
stagnation point of the pitot tube with 3 (see Fig. 1) 
and integrating between stations | and 3, Eq. (3) be- 
comes, with u; = 0, 


(4) 


Introducing the relation between the enthalpy 7 and 
the speed of sound ¢ 


[k/(k — 1)] (p/p) = — 1) 
Eq. (4) becomes 


= 2(43 11) 


(5) 


[2/(k — 1) Jer*[(es/a1)? — 1] (6) 


Substituting in this equation the temperature ratio 
for the square of the ratio of speeds of sound in accord- 
ance with Eq. (2), the following relation between the 
flight Mach Number J/; = u,/c, and the temperature 
ratio is obtained. 


M, = {[2/(R — 1)] — (7) 


By making use of the adiabatic relation between the 
temperatures and the pressures, the subsequent rela- 
tions between the flight Mach Number J, and the 
pressure ratios p3/p, and Ap/p, = (ps — pi)/p, follow. 


The dependence of the temperature ratio and the pres- 


uy" = 
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sure ratio on the flight Mach Number is then 
T3/T; =1+ rt 1)/2] (9) 


(B) Supersonic Air Speed 


In the supersonic case the equations developed for 
the subsonic case are not applicable but must be some- 
what modified. 

Since the nose angle of a pitot tube will always be 
above the critical angle because of the finite diameter 
of the pitot orifice, a bow wave will always precede the 
pitot tube when it travels with supersonic speed 
through the air. (See Fig. 1.) This bow wave is a 
curved detached shock wave that at some distance be- 
hind the pitot tube blends into the Mach cone. Fig. 2 
is an interesting schlieren photograph of a so-called 
\-shock in a wind-tunnel test. It clearly shows the 
pitot tube projecting into the A-shock, and the bow 
wave ahead of the nose of the tube. This photograph 
was taken at a Mach Number of 1.225 (Fig. 4 of refer- 
ence 10). 

The problem now is to find first the flow conditions 
through the shock wave and then the slowing up of the 
air stream to zero velocity at the stagnation point at 
the nose of the pitot tube. 

The only part of the bow wave that is of interest 
here is the center region through which the stream tube 
passes, which leads to the pitot orifice. This part of 
the shock wave will be curved, but the radius of curva- 
ture for this small area of the bow wave will be ex- 
tremely large, at least for moderately high Mach 
Numbers, so that it is permissible to assume here the 
radius of curvature to be infinite. Since this small 
area in the center of the bow wave will be perpendicular 
to the stream tube, the problem can be treated as if a 
plane shock wave were ahead of the pitot opening. 

From the general theory of plane compression shocks, 
it follows directly that the flow in the straight 
stream tube between the bow wave and the pitot orifice 
will be subsonic over its entire length (Prandtl’s rela- 
tion #2. = c, in which 2 designates the state behind 
the shock wave and c,, the critical speed of sound). 
Furthermore, no deflection of the velocity vector will 
occur because of the flow through the shock wave, so 
that the assumptions (d) (one-dimensional flow) and (f) 
(irrotational flow), which were made for the subsonic 
case, are still valid. These last assertions are valid, 
however, only for the small stream tube in the center 
of the shock wave. 

While the air cannot be considered an ideal gas with 
no viscosity in the bow wave proper, this assumption 
can be retained for the air before and behind the bow 
wave. The same holds true with regard to the adi- 
abatic flow before and behind the bow wave. 

The assumptions (a) (the air being a continuous 
medium), (e) (coordinate system moving with the ait- 
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MACHMETERS FOR HIGH-S 
TABLE 1 

Constants of Eqs. (15) and (16) for Various Values of k = cp/ce 
k A B D F 

1.4 166.9216 7 5/2 1.2876 
1.3 1628.218 26/3 10/3 1.2177 
1.2 285311.7 12 5 1. 1466 
1.1 2.852826 xX 103% 22 10 1.0741 


plane and steady flow condition), and (g) (no gravity), 
can also be retained. 

Condition (i), (constant entropy) cannot hold true 
through the shock wave. Also, the entropy will not 
be constant behind the curved shock wave, but will 
vary from point to point. However, the flow ahead of 
the bow wave and in the small stream tube behind the 
center part of the bow wave can be considered as being 
isentropic, but with different values of entropy. 

The temperature ratio between the stagnation tem- 
perature 73 and the free-stream temperature 7; for 
supersonic flow is identical with that for subsonic flow 
as given in Eq. (9). This follows immediately from the 
integration of the energy equation, Eq. (3). 


(1/2)us? + = (1/2)t2? + te = (11) 


In this equation suffixes 1 and 2 denote arbitrary states 
along a stream line, and the suffix 3 denotes the state of 
rest. The energy equation holds true whether states 
land 2 are on the same side of a shock or on opposite 
sides. The enthalpy of the state of rest, 7; is, there- 
fore, the same in both cases. Since for an ideal gas the 
enthalpy is 7 = c,7, and c, is a constant, it follows that 
the stagnation temperature 7; is independent of any 
shock condition. * 

The pressure ratio through the plane shock can be de- 
termined by applying to the flow the laws of conserva- 
tion of mass, of conservation of momentum, and of con- 
servation of energy. Designating the conditions behind 
the shock with the suffix 2, the pressure ratio through 
the shock becomes 


pe ( 2k ) 

k+1 (12) 
This relation was first given by Prandtl. The second 


law of thermodynamics—which states that the entropy 
ofa closed system can only remain constant, or increase, 
but cannot decrease—limits this equation to pressure 
imcreases through the shock—i.e., only compression 
shocks are possible. 
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After passing through the shock wave, the air is 
slowed down to zero velocity at the stagnation point 
at the orifice of the pitot tube, thereby being com- 
pressed adiabatically. The equations derived for the 
subsonic case can be applied to this compression, ex- 
cept that the conditions of state 2 behind the shock 
wave have to be substituted for the free-stream condi- 
tions of state 1. The adiabatic compression behind the 
shock wave is then 


= {1 + [(k — 1)/2] (Me) (13) 


in which M, represents the Mach Number behind the 
shock wave, which is smaller than 1.0. 

Combining the two pressure ratios, Eqs. (12) and 
(13), Lord Rayleigh, in 1910, obtained the following 
equation for the ratio of the stagnation pressure to the 
static pressure :!! 


ps _ ) (= 
k+1 


pi 
1 k/(k—1) 
(14) 


+ 1)? - ok 


4k a 


To extract the flight Mach Number /, in a simple 
exact form from this somewhat cumbersome equation 
seems a hopeless task. However, the equation can be 
simplified somewhat by introducing the magnitude of 
the ratio of the specific heats k. Then it can be written 
in the form 


= — (1/M,*)}? (15) 


With increasing Mach Number this equation tends to 


the limiting equation 
(ps/P1) = (16) 


The constants A, B, D, and F, which are dependent on 


k, are listed in Table 1. For air, with k = 1.400, 
Eqs. (15) and (16) read 

= 166.9216 M,2/[7 — (17) 
and 


(Ps/Pi) = 1.2876 Mi? (18) 


For a Mach Number of 4.0 in air, the difference be- 
tween the two equations is approximately 2 per cent. 
(See Table 2.) 


TABLE 2 


Pressure Ratio: Stagnation Pressure to Static Pressure p;/p; for Various Values of k = cp/c, and Different Flight Mach Numbers 1; 
and the Differences Between the Values p;/p; in Percentages Obtained by Using the Limiting Eq. (16) and the Exact Eq. (15). 


k = 1.4 ~ k= 13 ———-k = ———-k = 
M ps/ Diff. % Diff. % Diff. % pr Diff. % 
1.0 1.893 — 32.0 1.832 —33.5 1.772 —35.3 1.710 —37.2 
1.5 3.413 —15.1 3.265 —16.1 3.116 —17.2 2.964 —18.5 
2.0 5.640 —8.7 5.370 -—9.3 5.096 —10.0 4.817 —10.8 
3.0 12.06 -—3.9 11.44 —4.2 10.81 —4.5 10.17 —-4.9 
4.0 21.07 —2.2 19.96 - —2.4 18.83 —2.6 17.68 -—2.8 
9.0 32.65 —1.4 30.92 -—1.5 29.15 —1.7 27.35 —-1.8 
6.0 46.82 —1.0 44.31 -1.1 41.76 —-1.2 39.16 —-1.3 


4 
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Fic. 3. Diagram of ratios: stagnation pressure to static 
pressure p;/p, and stagnation temperature to static temperature 
T;/T, vs. flight Mach Number /, for various values of k = 
Cy | 


T; 
T, 


Zahm derived an expression somewhat similar to 
Eq. (15), which, however, is not quite so simple as the 
above.” 

For high Mach Numbers these equations are not 
valid anymore, because appreciable heat exchange will 
occur, and condition (h) (adiabatic flow), under which 
the above equations were derived, no longer holds true. 
They are also not valid for very high altitudes, for there 
the air can be considered no more as a continuous 
medium, because of the long molecular paths due to the 
low density. 

The curves for the temperature ratio 73/7, and the 
pressure ratio p;/p; versus the free stream, or flight 
Mach Number, are shown in Fig. 3. Since, in the near 
future, there may be some interest in Machmeters for 
wind tunnels that are operated with a different gas than 
air, families of curves for different values of the ratio 
of the specific heats k have been plotted. These curves 
show that a standard Machmeter, calibrated for air, 
can be recalibrated without too much difficulty for 
some Other gas. For instance, the N.A.C.A. has investi- 
gated the possibility of using Freon-12 for aerodynamic 
testing. This gas has a value of k = 1.125, its density 
is approximately four times that of air, and its velocity 
of sound is less than one-half that of air.'” 


Although the equations of the pressure ratio for the 
subsonic and the supersonic cases are different, the pres- 
sure ratio curves are continuous at Mach Number J, = 
1.0, since the first derivatives of both equations, which 
indicate the slopes of the curves, are identical at this 


point. However, the second derivatives, of which the 
radii of curvature are functions, are different for both 
equations at that particular point. 

Fig. 4 shows the actual calibration curves of Mach- 
meters operated by the differential pressure Ap = 
ps — pi fork = 1.40. A family of curves of Ap versus 
flight Mach Number M, for altitudes h from 0 to 
100,000 ft. has been plotted. The barometric pressures 
for the various altitudes have been taken to be the 
pressures of the U.S. Standard Atmosphere.'* It is 
at once apparent that such a Machmeter must have an 
altitude correction or, expressed more precisely, a baro- 
metric pressure correction. The mechanism of a Mach- 
meter consists, therefore, of a lever arrangement actu- 
ated by a diaphragm that responds to the differential 
pressure; with changes in altitude, this lever arrange- 
ment is modified by a second diaphragm that responds 
to barometric pressure changes. 


(V) THe Sussonic MACHMETER 


The subsonic Machmeter had already been developed 
during the war for the Army. It has been used ex- 
tensively, both in test work and in flight operations. It 
is calibrated in accordance with Eq. (8) for a range of 
flight Mach Numbers 0.3 to 1.0 and a range of altitudes 
from 0 to 50,000 ft. (see Fig. 4) with an accuracy of 
+().01 Mach Number over the entire range of calibra- 
tion. The temperature range of operation is from 
—55° to +70°C. with an accuracy of +0.02 at these 
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Fic. 4. Diagram of differential pressure Ap = p3 — p, vs. flight 
Mach Number MV, for various altitudes h and k = 1.40. 
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extreme temperatures. More accurate calibration is 
possible but requires a great deal more calibration 
time than the standard instrument. 

The differential pressures operating this instrument 
range from 0.221 to 26.718 in. Hg (0°C.), the first value 
corresponding to M1, = 0.3 at A = 50,000 ft., while the 
last value corresponds to M, = 1.0 at sea level. 

A schematic diagram of the operation of the subsonic 
Machmeter is shown in Fig. 5. The stagnation pressure 
or pitot pressure p; is led into the open diaphragm, while 
the static barometric pressure /; is led into the inside 
of the airtight case of the instrument. The resulting 
movement of the open diaphragm, which is a function 
of the differential pressure Ap, by means of a lever arm 
rotates the first rocking shaft, which in turn rotates the 
sector rocking shaft through a pair of levers. From 
this sector rocking shaft the movement is transmitted 
to the handstaff by means of a sector and pinion. The 
altitude correction is accomplished with an evacuated 
altitude diaphragm, which responds to the static pres- 
sure ~;. With increase in altitude—i.e., decrease in 
static pressure—this diaphragm expands and pushes 
the first rocking shaft to the left. This movement 
causes a shortening of the effective lever arm of the 
sector rocking shaft lever, thereby changing the amount 
of rotation caused by the open diaphragm. 

The temperature compensation of this instrument is 
accomplished with a bimetallic bracket mounted on the 
altitude diaphragm, as is shown in the diagram. This 
bracket preloads the diaphragm by means of two push 
pins, which engage the center piece, thereby changing 
the expansion of the diaphragm in such a manner as 
to offset the temperature error of the instrument. 

The actual instrument is not so simple as this sche- 
matic diagram may indicate. The open diaphragm 
is actually restrained by a leaf spring; screws restrain 
the deflection of this spring so that with increase of 


t 


Fic. 5. Schematic diagram of the subsonic Machmeter. 


Schematic diagram of the transonic Machmeter. 


Fic. 6. 


Mach Number the pointer movement is uniform. The 
altitude diaphragm must be made with extreme accu- 
racy so as to meet a certain expansion characteristic 
with change of barometric pressure. Furthermore, the 
two levers connecting the two rocking shafts are not 
perpendicular to their respective shafts, as shown in the 
schematic diagram, but one of them is in an oblique 
position. In spite of a number of adjustments, which 
are also not shown in the diagram, the calibration of 
this instrument is rather delicate and tedious, since any 
adjustment made on one part of the mechanism has 
immediate effect on the other part. 


(VI) THe TRANSONIC MACHMETER 


The transonic Machmeter, which is rather a recent 
development for the present high-speed research, is 
calibrated in accordance with Eqs. (8) and (17) for the 
subsonic and supersonic speeds, respectively. Its 
range is for flight Mach Numbers 0.5 to 1.5 and for 
altitudes from 0 to 50,000 ft. The first instruments 
had an accuracy of 0.01 Mach Number for the range 
of 0.65 to 1.35; the two extremes of the range were 
within +0.015. It is expected that future instruments 
will have the accuracy of +0.01 throughout. The 
temperature range of operation is the same as that of 
the subsonic Machmeter. 

The differential pressures operating the instrument 
are from 0.640 to 72.208 in. Hg (0°C.), the first value 
corresponding to M, = 0.5 at h = 50,000 ft., while 
the last value corresponds to MM, = 1.5 at sea level. 

The first intention was to use, for this transonic 
Machmeter, a modified mechanism of the subsonic in- 


strument. Preliminary tests with restraining of the dia- 
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phragm, so that a uniform scale on the dial could be 
obtained, did not look too promising because of the 
high maximum pressure of 72.2 in. Hg (2.7 Atm.). It 
was then decided to build a mechanism with a dia- 
phragm, the deflection of which is linear with pressure 
increase, and which would not be restrained at all. 
The diaphragm of the 100 in. Hg Manifold Pressure 
Gage, which has good general characteristics and which 
was available, was chosen for this purpose. The 
straightening out of the strongly curved differential 
pressure curve versus Mach Number (see Fig. 4) was 
to be done by a pair of oblique levers. It was also felt, 
that, although this new mechanism would be some- 
what more complicated, the calibration would be sim- 
pler and would require less time. 

A schematic diagram of the first transonic Mach- 
meter is shown in Fig. 6. The positions of the levers 
shown with full lines represent the condition of zero 
differential pressure and sea-level altitude; those indi- 
cated by dash-dot-dash lines show extreme positions 
for differential pressure and altitude. 

The operation of the instrument is as follows: 

As in the case of the subsonic Machmeter, the stag- 
nation or pitot pressure /; is led into the open diaphragm 
E, while the static pressure /; is led into the inside of the 
airtight instrument case. The movement of the open 
diaphragm is transmitted by means of the lever F, 
which is mounted on rocking shaft A, to rocking shaft 
B. The rotation of this rocking shaft is essentially pro- 
portional to the differential pressure. Then a pair of 
oblique levers J and H transmits the rotation of rock- 
ing shaft B to rocking shaft C; the obliqueness of the 
levers is so computed that the rotation of rocking shaft 
C is approximately proportional to the Mach Number. 
The pointer of the instrument is rotated by means of a 
sector on rocking shaft C and a pinion on the handstaff 


D. 


Fic. 7. Front view of the transonic Machmeter. 
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The altitude correction is attained in the following 
manner: 

Each calibration curve of the family of curves in 
Fig. 4 is spaced from the abscissa in the ratio of the 
barometric pressure at any altitude to the sea-level 
pressure, p,/ps,;. This means that, for instance, if an 
aircraft is flying at a Mach Number of 1.2 at an alti- 
tude of 50,000 ft., the deflection of the air-speed di- 
aphragm will correspond to a differential pressure of 
only (3.437/29.921)42.114 = 4.838 in. Hg; without 
an altitude correction the pointer would, therefore, 
indicate a Mach Number of about 0.46. In order to 
read the correct Mach Number of 1.2, the rotation of 
rocking shaft B has to be speeded up in the proper ratio. 
By moving rocking shaft A closer to rocking shaft B, 
the effective length of the lever arm G is shortened. 
For a certain rotation of rocking shaft A, the rotation 
of rocking shaft B is now increased. This axial displace- 
ment of rocking shaft A is accomplished by the altitude 
diaphragm J, which expands with decrease of static 
pressure. 

The actual instrument is not so simple as the sche- 
matic diagram. A set of four cells is used for the open 
diaphragm, which has a straight-line movement with 
increase in pressure. The total deflection of the dia- 
phragm for the full range of 1.5 Mach Numbers is 
0.260 in., while the deflection for the lowest calibration 
point of Mach Number 0.5 is only 0.0023 in. 
small movement deflects the pointer approximately 
140° from the zero point of the mechanism to the gradu- 
ation 0.5. On the finished instrument, however, the 
pointer is restrained by a stop at the zero graduation, 
which is only about 25° away from the 0.5 graduation. 
Fig. 7 shows the dial arrangement of the first transonic 
Machmeter. 

The specification called for a dial with nearly uni- 
form spacing of the graduations. This meant that the 
oblique lever arrangement of.levers J and H would have 
to furnish, for the low range, a considerable amplifica- 
tion of the diaphragm movement, since the curve, 
differential pressure vs. Mach Number (Fig. 4), 
has a small slope in that part of the range; with in- 
crease of Mach Number, this amplification had to 
decrease rapidly. The lever arrangement, which was 
determined theoretically, seemed to have, for the low 
range of Mach Numbers, unfavorable angles between 
the levers as far as the friction was concerned. A mock- 
up, built subsequently, showed the lever relations to be 
satisfactory, and the dial spacing over the entire range 
was as nearly uniform as that of the supersonic region 
of the final dial shown in Fig. 7. When the mechanism 
itself was built and tested, however, it was found that 
the lever arrangement was a little too close to the criti- 
cal friction angle. This might have caused jamming of 
the levers under adverse conditions. A slight change 
in the angular relation between the levers J and /I 
removed this potentially dangerous condition; but it 
also somewhat spoiled the appearance of the dial in 
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Fic. 8. Schematic diagram of the altitude compensation of the 
transonic Machmeter. 


regard to uniformity at the beginning of the scale, as 
can be seen in the photograph. 

An interesting feature of the instrument is the 
mechanism for the altitude correction. It is shown 
schematically in Fig. 8. 

Rocking shaft A, with its lever F and counterweight, 
is jewel-mounted in a balanced magnesium frame K. 
The frame is supported by two highly polished and 
hardened slide rods L and M of 0.031-in. diameter. 
The ends of these rods, which are clamped into the 
frame A, have 0.012-in. diameter pivots carrying the 
rocking shaft A. The rods slide in jewel bearings. In 
spite of the considerable complication, rocking shaft A 
is mounted in separate bearings of small size in order to 
reduce the friction of the mechanism to a minimum. 
The axial movement of this assembly is accomplished by 
a conical hairspring NV, which is mounted on a collet 
clamped on the right end of the slide rod L. This 
conical hairspring has a dual purpose. First, it is in- 
stalled in such a manner that it will always try to pull 
the rocking shaft assembly to the left. The altitude 
diaphragm controls the amount of axial displacement 
the spring can effect. Second, the hairspring turns the 
sliding rods while an axial displacement takes place. 
This rotation, which is approximately 35° for an axial 
displacement of slightly more than | in., reduces the 
friction in the comparatively large jewels in which the 
sliding rods are mounted. 

The altitude diaphragm J, which affects the altitude 
compensation, consists of two evacuated cells. Its 
expansion for an altitude change from sea level to 
50,000 ft. is 0.320 in. This movement is transmitted 
through a crank mechanism to rocking shaft O on which 
a sector P with a groove at its periphery is mounted. 
Connection between this part of the mechanism and 
the rocking shaft A is made by an extremely thin beryl- 
lium copper ribbon Q, which is attached on one end to 
the hub of the conical hairspring N and at the other end 
to the sector P after it is guided through its groove. 
The expansion of the altitude diaphragm is matched to 
the crank mechanism motion in such a manner that the 


resulting movement displaces rocking shaft A in the 
proper amount so as to change the effective lever arm 
of lever G for the required altitude compensation. 

Fig. 9 shows a photograph of the delicate assembly of 
rocking shafts A and B on a special mounting plate R. 
The position of the various components corresponds 
approximately to a Mach Number of 1.5 at an altitude 
of 10,000 ft. The magnesium frame A, the rocking 
shaft A with its lever F and counterweight, the sliding 
rods L and M, and the conical hairspring NV can easily 
be recognized. The beryllium copper ribbon Q, one 
end of which is attached to the collet of the hairspring 
while the other is free, can also be noted. The little 
brass collet W near the jewel bearing of the slide rod 
M has been installed only in order to hold the assembly 
in this position while being photographed. If it were 
removed, rocking shaft A would be pulled by the hair- 
spring all the way to the left beyond the 50,000-ft. 
position. 

The engagement between levers F and G is clear on 
this photograph. Lever H/ shows a slight bend near the 
rocking shaft B. This bend represents the slight change 
in the lever arrangement between levers /7 and /, which 
improved the friction conditions considerably. 

It is important that the slide rod S, from which the 
lever F picks up the diaphragm movement, is perfectly 
parallel to the rocking shaft A. For this reason, the 
slide rod S is mounted on a separate frame 7. The 
pivots about which this frame turns are inserted in the 
same mounting plate R, in which the pivots and jewels 
for rocking shafts A and B are mounted. The frame 
rocks only over a small angle, so that the lateral dis- 
placement of slide rod S is small and can be compen- 
sated by the calibration of the instrument. The frame 
is carefully balanced by the counterweight U. The 
locations of the pivots are indicated by the letters V. 
The differential pressure diaphragm F engages the slide 
rod S by means of another rod that is resting on it. 


Fic. 9. Slide mechanism of the transonic Machmeter. 
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Fic. 10. Side view of the transonic Machmeter with the instru- 
ment case partially cut away. 


The whole slide mechanism, consisting of rocking 
shafts A and B and frame 7, as seen in Fig. 9, is placed 
into the mechanism; the proper line-up with the rest 
of the instrument is secured by two locating pins. This 
slide mechanism can be seen in place in Fig. 10, which 
shows a side view of the instrument with most of the 
case cut away. 

The mechanism frame consists of a front plate, a 
back plate, and three long posts connecting them. 
On two of the posts the two sets of diaphragms are 
mounted, while the third carries the slide mechanism. 
The sector rocking shaft O of the altitude mechanism 
is mounted in jewels in a special bracket, which can be 
seen on the left side of the photograph. The odd-shaped 
part covering most of the sector rocking shaft is the 
counterweight for the sector and the lever arm. The 
rocking shaft C is mounted between the front plate and 
a post, which is inserted into one of the larger posts. 
The oblique lever arm / of this rocking shaft can be 
seen in engagement with the lever arm H of rocking 
shaft B on the right side of the mechanism. 

Both diaphragms are weight compensated, so that 
errors due to acceleration are negligible. Of course, 
all movable parts are carefully balanced for the same 
reason. The temperature compensation of this instru- 
ment is done on the altitude mechanism. 

The size of the transonic Machmeter is the con- 
ventional 3!/s-in. diameter, with a case length of 6 in. 
The connections in the back are identical to those of a 
conventional air-speed indicator. 


(VII) THe Macu NuMBER Limit SwIitcH 


The two Machmeters described above are used 
chiefly in piloted aircraft where a continuous indica- 
tion of the flight Mach Number is required. However, 
for a pilotless missile, such as the buzz bomb, a con- 
tinuous indication is not necessary. Here the only re- 
quirement is to prevent the missile from flying faster 
than its critical Mach Number. The so-called Mach 
Number Limit Switch, a special version of a Mach- 
meter, performs this task. This instrument opens and 
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closes electrical circuits within specific limits that can 
be set from the outside of its case. So far, there is only 
a subsonic model of the Mach Number Limit Switch, 
which is operating in accordance with Eq. (8). If the 
demand arises, this range can be extended into the 
supersonic domain. 

The instrument contains, as the other Machmeters 
do, an open differential pressure diaphragm and an 
altitude diaphragm. Both diaphragms rotate levers to 
which electrical contacts are attached. At the proper 
Mach Number settings, these contacts are opened or 
closed, thereby controlling a servomotor that changes 
the rate of the fuel supply to the engine. 


(VIII) THe Macu AIRSPEED INDICATOR 


A paper on Machmeters would not be complete with- 
out a short discussion of its sister instrument, the 
Maximum Allowable Airspeed or Mach Airspeed Indi- 
cator. This instrument, which was developed for opera- 
tional use, was used extensively during the war on 
fast fighter aircraft and is now being increasingly 
adopted in fast-flying commercial transport aircraft. 

The Mach Airspeed Indicator is a combination in- 
strument consisting of two entirely different mecha- 
nisms, the pointers of which indicate on the same dial. 
One of these instruments is a conventional air-speed 
indicator with approximately one revolution of the 
pointer for the full range. The mechanism of the other 
instrument, which operates a warning pointer, is actu- 
ally an altimeter that is calibrated in maximum allow- 
able air speed for a certain Mach Number. 


Before installation in a particular aircraft, this Mach 
mechanism, as it is usually called, is set by means of a 
screw driver to the highest flight Mach Number at 
which this type of aircraft can safely be flown. This 
Mach Number appears in a small window of the dial. 
The indicated air speed corresponding to the safe Mach 
Number is then indicated by a warning pointer, which is 
somewhat wider than the air-speed pointer; the warn- 
ing pointer is painted either red or in red and yellow 
stripes. As the aircraft ascends, the warning pointer 
indicates lower and lower air speeds, because the indi- 
cated air speed for a given Mach Number diminishes 
with altitude, since the value of the speed of sound de- 
creases with increase in altitude. 

The mathematical relation between the indicated 
air speed V; and the altitude / for a given flight Mach 
Number /, for subsonic speeds follows from the sub- 
sequent considerations: 

The differential pressure AD = p; — Pp; for any indi- 
cated air speed is a constant, regardless of altitude. 

(19) 


Ap, = Apsz, for V; = const. 


The subscripts designate any altitude aad sea level. 
By substituting Eq. (10) with the appropriate suf- 
fixes, Eq. (19) reads, with p, and ps, as the static 
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Fic. 11. Diagram of altitude h vs. indicated air speed V; for 
various flight Mach Numbers M, and k = 1.40. 


pressures and M, and Ms, as the free-stream Mach 
Numbers at any altitude and at sea level, 


+ (2) i} (20) 


Under standard sea-level conditions, the true air 
speed equals the indicated air speed, so that 


Vsz = Vi (21) 
Substituting the Mach Number at sea level 
Msz = (V/c)st = (22) 


in Eq. (20) and solving for the indicated air speed, the 
following relation is obtained : 


k/(k—1) (k—1)/k % 


This value represents the Maximum Allowable Indi- 
cated Airspeed for any one flight Mach Number as a 
function of the static pressure p, at any altitude Fig. 
11 is a graph showing the relation between indicated 
air speed V;, and altitude A for various flight Mach Num- 
bers M,. 

The calibration of the air-speed mechanism is in 
accordance with the well-known relation developed 
more than a hundred years ago by de Saint-Venant and 


Wantzel in 1839. 
Qk Pst 

ti — lpsz Psi (24) 


This equation is obtained by setting the right-hand 
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side of Eq. (20) equal to the differential pressure Ap and 
introducing the indicated air Speed V; by means of 
Eqs. (22) and (2). 

Fig. 12 shows a schematic view of both mechanisms 
of the Mach Airspeed Indicator. The stagnation or 
pitot pressure ~; is led into the open diaphragm, while 
the static pressure is led into the airtight instrument 
case. The expansion of the open diaphragm is, there- 
fore, a function of the pressure differential Ap = p3 — 
pi. 

The air-speed mechanism is operated by the open 
diaphragm A, which rotates the rocking shaft B by 
means of a crank mechanism. Rocking shaft B in 
turn, by means of a sector and pinion, rotates the 
handstaff C, which carries the air-speed pointer D. 

The Mach mechanism is operated by a closed evacu- 
ated diaphragm £, the expansion of which is a function 
of the static pressure, which rotates rocking shaft F by 
means of a sliding lever. The rotation is transmitted 
from there to rocking shaft G through a pair of right 
angle levers and then, by means of a sector and pinion, 
to the handstaff H, which carries the warning pointer 
I. The change in calibration for the various flight 
Mach Numbers is accomplished through rotating the 
shaft J, which carries a cam and which is connected 
to the subdial K by gears. Rotation of the cam will 
effect an axial displacement of rocking shaft F, thereby 
changing the effective length of the lever arm on rock- 
ing shaft G. 

As an additional feature, the instrument has an ad- 
justable stop L, which may limit the movement of the 
warning pointer. This feature is of importance at low 


Fic. 12. Schematic diagram of the Mach Airspeed Indicator.. 
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Fic. 13. Schematic diagram of two supersonic static pressure 
probes. 


altitudes, where the placard speed of the aircraft may 
be lower than the maximum allowable air speed with 
respect to compressibility dangers. 

The Mach Airspeed Indicator is used by the pilot in 
the same manner as the conventional air-speed indi- 
cator. Simultaneously, in addition to the air-speed 
indication, he can see from the angle formed between 
the two pointers how near to the dangerous air speed 
his flying speed is; as long as the two pointers do not 
cross, he is flying safely. 

The Germans had a different type of Mach air-speed 
indicator. This instrument was a conventional air- 
speed indicator for low speeds and a Machmeter for 
high speeds. They probably used it in their rocket- 
propelled Me 163 fighters. 


(IX) LocaTION OF PRESSURE VENTS 


The location of the pressure vents for transonic and 
supersonic air speeds represents a number of problems. 
For the measurement of the stagnation pressure, pitot 
tubes of shapes similar to those used for subsonic flow 
will most likely be usable, as long as they are placed far 
enough from any possible interference effects resulting 
from shock waves. 

The proper location of the static vent is a more 
difficult problem. The use of conventional pitot tubes 
is not feasible because a bow wave forms in front of the 
nose of the tube, and the pressure behind that bow wave 
is higher than that in the undisturbed atmosphere; 
locating the static vent holes near the nose of the air- 
craft at the side of the fuselage, as is done on the P-80 jet 
airplane, suffers from the same evil. Pressure probes 


of the type used for wind-tunnel tests may not readily 
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be used for aircraft because of their extremely small 
diameter of about 1 mm.!° 

It is conceivable that the atmospheric pressure can 
be obtained by using a static probe that has the vent 
hole on the straight side of a knife edge (Fig. 13a). In 
this arrangement a shock wave would be created at the 
side that is beveled, while a Mach line, which repre- 
sents only a very weak shock, would occur on the 
straight side of the knife edge. If such a static pressure 
sonde could be aligned perfectly with the oncoming 
air stream, the ambient atmospheric pressure could be 
obtained. 

It seems feasible to design a static pressure probe for 
supersonic flows by bending the knife edge into a cylin- 
der (Fig. 13b); such a probe would look somewhat like 
a Venturi tube, except that a number of pressure vents 
would be located evenly spaced on the outside. If 
properly lined up with the air stream, a shock wave 
would form, which would be drawn into the tube open- 
ing; on the outside, only a Mach cone would result, 
the disturbance of which is negligible, so that the am- 
bient atmospheric pressure could be measured. The 
rim angle would have to be extremely small and the 
inside diameter of the probe, comparatively large. The 
static probe could be kept in line with the air stream by 
servomotors; the differential pressure between the four 
static vents, which would exist when the probe was not 
in proper alignment, could control the servomotors. 

The location of the static pressure vents somewhere 
on the fuselage is another possibility for static pressure 
determination. Locations on the surface of aerody- 
namic bodies, where the surface pressure is identical with 
the atmospheric pressure, also exist in supersonic flows. 
For instance, wind-tunnel measurements of the pressure 
distribution over the A-4 rocket (V-2) show that for 
free-stream Mach Numbers from 1.56 to 3.24 there is 
such a location about 10 m. aft of the nose. Calibra- 
tion tests up to a Mach Number of 2.0 on actual A-4 
rockets in White Sands, N.M., show excellent agree- 
ment between the measured values of atmospheric 
pressure during the rocket ascent and values obtained 
from balloon sondes at about the same time. Only 
in the range near Mach Number 1.0 are the values taken 
by the A-4 rocket somewhat higher than those of the 
balloon sonde.” Such locations where the surface 
pressure is identical with the atmospheric pressure can 
probably also be found on aircraft flying at supersonic 
speeds. 


(X) MacCHME’ER BASED ON TEMPERATURE 
MEASUREMENT 


The use of temperature as the parameter to operate 
a Machmeter is possible in principle. The mechanism 
could be made similar to that of the pressure-operated 
instrument. The diaphragms would now be connected 
to temperature bulbs, which would be installed outside 
the aircraft: one, for measuring the stagnation tem- 
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perature and the other, for measuring the free-stream 
temperature. Liquid in the temperature bulbs and the 
connecting capillaries would be the medium of opera- 
tion. 

From the point of view of instrument design, such a 
scheme seems perfectly feasible. However, the one 
difficulty lies in the measuring of the proper values of 
the stagnation temperature 7; and the atmospheric 
temperature 7, at transonic and supersonic speeds, 
or at least measuring values that can be corrected in a 
simple manner. The measurement of these tempera- 
tures, however, represents a difficult task because of 
radiation and convection effects. These are functions 
of the altitude and the Mach Number, the latter being 
the parameter to be obtained from the instrument. In 
a report by Wood," a graph showing equilibrium sur- 
face temperatures as a function of free-stream Mach 
Number for various altitudes indicates that the correc- 
tions that would have to be applied to the actually 
measured temperature values would be extensive. 
These considerations seem to rule out such Machmeters. 


(XI) CoNncLusIoNns 


Machmeters will acquire more and more importance 
for the steadily increasing flying speeds. Greater range 
in Mach Number, as well as in altitude, will be needed 
in the near future. Of these two problems, the increase 
in altitude range represents the greater difficulty. One 
glance at Fig. 4 or Table 3 shows the small differential 
pressures available at 100,000 ft., for example, as com- 
pared with those at sea level. 

It is likely, that for high-speed aircraft the conven- 
tional air-speed indicator will be superseded by the 
Machmeter. Perhaps a combination instrument, 
showing indicated air speed for the low-speed range 
and Mach Number for the high-speed range, might 
eventually be the best solution for a flight speed instru- 
ment. 

The location of the static vent is one problem in 
connection with Machmeters which still seems to re- 
quire considerable experimental research, particularly 
at high Mach Numbers. 

From the experience gained in developing several 
Machmeters, it seems certain that, as demands for 
greater ranges in Mach Number and altitude occur, new 
Machmeters can be developed to keep step with those 
requirements. 


TABLE 3 


Differential Pressure Ap = ps — p, (In. Hg 0°C.) for Various 
Mach Numbers /, and Altitudes 4 (Ft.) 


M, h =0 50,000 100,000 150,000 

pi = 29.9212 3.4370 0.3225 0.0444 
0.5 5.57 0.640 0.0601 0.0083 
1.0 26.718 3.069 0. 2880 -0.0396 
1.5 72.208 8.294 0.7783 0.1071 
2.0 138.847 15.949 1.4965 0.2060 
2.5 225.191 25. 867 2.4272 0.3342 
3.0 330.957 38.017 3.5672 0.4911 
3.5 456, 059 52.387 4.9155 0.6767 
4.0 600.461 68.974 6.4720 0.8910 

k = 1.400. 
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An Improved Schlieren Apparatus’ 


FRANK W. BARRY? anp GILBERT M. EDELMAN? 
Massachusetts Institute of Technology 


SUMMARY 


This paper describes an improved schlieren apparatus that 
forms two separate images of the test section simultaneously. 
One image indicates vertical density gradients; the other, 
horizontal density gradients. The determination of transient 
density fields is more accurate with this modified schlieren 
system than with the ordinary schlieren system taking two 
successive pictures. At present, this apparatus is being used 
to study the interaction between shock waves and boundary 
layers. 


DISCUSSION 


SCHLIEREN APPARATUS is an optical instru- 


ment that is used in the study of supersonic flow 
phenomena. With this instrument photographs are 
taken in which the pattern and intensity of the light 
and dark areas are a measure of the density gradients 
in the field of view. 
In the usual construction of the schlieren equipment, 
a single straight knife edge is used, and the resulting 
picture shows density gradients normal to that knife 
edge. If density gradients in another direction are 
desired, the knife edge must be rotated and another 
picture taken. 
The time required to make such an adjustment, and 
hence the time lag between successive pictures, is not 
significant when the flow is steady. However, when 
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the flow is not steady and schlieren pictures corre- 
sponding to two different knife edge orientations are 
desired, such pictures should be taken simultaneously 
to obtain comparable pictures. 

A device has been developed at the Gas Turbine 
Laboratory of the Massachusetts Institute of Tech- 
nology for this purpose. Two pictures are taken 
simultaneously on a single negative. One of these 
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pictures shows density gradients in a vertical direction; 
the other shows density gradients in a horizontal direc- 
tion. In this device the light from a single line source 
is split by optical means into a horizontal and a vertical 
beam. Each of these beams passes one of two mutually 
perpendicular knife edges in the schlieren system and 
produces the two separate images on the viewing 
screen or film. Each beam is independent and the 
effect is such as would be produced by two separate 
schlieren systems. 

Fig. 1 is an illustration of the complete schlieren 
apparatus with the double knife-edge arrangement. 
The light source unit is shown on the left and the 
camera unit on the right. Light from source A passes 
through lens B to the light source converter, C. 

Fig. 2 illustrates the optics of the light source con- 
verter, C. The flat, converging light beam is split into 
two parts by an aluminized glass block. One of the 
two resulting beams emerges without change in either 
orientation or direction. The other beam emerges 
without change in direction, but in passing through 
the mirror system it is rotated through 90°. 

Referring to Fig. 1, the slits D are placed at the 
point where the two beams from C focus as an L- 
shaped image of the sources. The light that passes 
through the slits D is reflected by the plane mirror E 
to the parabolic mirror F, from which it emerges as a 
parallel beam composed of the light from both of the 
beams that emerged from C. This parallel beam passes 
through the test section G to a second parabolic mirror 
and a second plane mirror and forms a second L- 
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shaped image of the sources at the knife edges //. 
The knife edges H are similar in construction to the 
slits D, which are detailed in Fig. 3. 

From H, the two distinct beams pass to the image 
separator shown in Fig. 4. The image separator in- 
creases the distance between the beams and permits 
the formation of two separate images of the test section 
G on the film by the lens A. 

Fig. 5 is presented as an example of the photographs 
which have been taken with this apparatus at an ex- 
posure of about 5 microsec. The upper picture shows 
density gradients normal to the flow, such as exist in 
the boundary layers, and the lower picture shows den- 
sity gradients parallel to the flow direction. The 
flow direction is upwards. 


SHUTTER 


A Note on F light F lutter Testing 


ROBERT ROSENBAUM* anv ROBERT H. SCANLAN 
Republic Aviation Corporation 


SUMMARY 


An attempt has been made to correlate theoretically calculated 
decay effects of wing oscillations at subcritical flutter speeds with 
results obtained during actual flight flutter testing at these speeds. 
The theory discussed has been found to give reasonably good 
agreement with the limited number of test data obtained in flight 
from a single airplane. Conclusions of a general nature are as yet 
unwarranted, but at least a new point of view is presented which 
applies to interpretation of both test and theoretical results. 


INTRODUCTION 


in the progress of their work by a serious lack of 
practical confirmation of their calculated results. On 
the other hand, results of flight flutter testing have not 
been very well supported by theory; flight testing for 
flutter in aircraft must necessarily occur below critical 
flutter velocities where usual theory does not hold. 
Flight testing for a condition of flutter is a difficult and 
possibly dangerous procedure. The typical flight flutter 
testing that has been done has been aimed merely at 
predicting the actual flutter speed of the aircraft in 
question. The procedure is usually to measure and 
plot some parameter, such as the rate of decay of 
flutter-type oscillations versus true air speed; this plot 
is the usual flutter trend or g-V curve. As the rate of 
decay lessens, it is assumed that a truly harmonic or 
flutter condition is approached; extrapolation of the 
trend curve is then used to indicate approximately 
the velocity at which flutter will occur. Knowledge of 
the actual flutter speed is, of course, the primary engi- 
neering goal of the usual test, but this test has two 
shortcomings: first, the actual flutter speed may be 
approached rather abruptly and, hence, not safely; and 
second, even if approximated by extrapolation of test 
results, the critical speed represents only a single 
“test” item to be compared against a single calculated 
item. Usually insufficient attention is given to a com- 
parison of subcritical oscillatory test items with theory. 
The trend curve is thus, in a sense, used incorrectly in 
that the most engineering attention is concentrated 
not on the test points that comprise the curve but on 
the extrapolated flutter point that is off the curve. 
When, during tests, damped subcritical oscillations of 
flutter type occur, it would be of extreme practical 
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value if they could be used not only in the usual test 
trend curve (to be extrapolated later) but also in direct 
comparison with predictions from a trend curve given 
by some theory. Thus, more attention would be given 
to the best items of the trend curve—namely, those 
points on it—and perhaps more validity given also to 
its least accurate item—the extrapolated flutter veloc- 
ity. It is entirely conceivable, for example, that 
extremely small rates of decay can occur in subcritical 
flutter oscillations which nevertheless occur at veloci- 
ties far below the critical flutter velocity. The question 
is: What should be considered a small rate of decay? 
A single test oscillation showing an apparently small 
rate of decay at a given air speed might, in itself, be 
alarming unless placed in its proper position in a trend 
curve. But hesitancy to complete the trend curve by 
flying to test velocities beyond this point can be readily 
understood unless some theory exists to point out the 
noncritical nature of such a condition. 

The aim of this note is to present a simple modifica- 
tion of standard flutter theory, which suggests a possible 
means of prediction of preflutter decay rates over a con- 
siderable range (hence, also a safe range) of air speeds. 
This theory is supported by some test results that are 
in reasonably good agreement with it. The ultimate 
validity of this theory rests on support by many more 
such test results. 


THEORY! 


The typical flutter stability determinant, evaluated 
for a given flutter parameter 1/k and set equal to zero, 
results in a polynomial, the complex roots of which 
have a form like the following: 


2 


where g is the damping coefficient required for flutter 
at the frequency w and w, is (as taken in this case) the 
“torsional’’ natural frequency of the system considered. 


Thus, 
X = (w/w)? 
and the critical flutter velocity is given by 
= bywa(1/k)/WX 


or 


\ (1/k 
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A NOTE ON FLIGHT FLUTTER TESTING 


where }, is some reference semichord (in feet) and 1/k 
is the dimensionless flutter parameter v/b,w. Also, g 
is calculated directly from g = Y/X. 

Using a given set of 1/k values, g is plotted versus 
/Wq to present the flutter characteristics of the airplane 
in question in the particular flutter mode investi- 
gated. 

It is assumed in obtaining the above g-V curve that 
motion of the system is neither damped nor diverging, 
but simple harmonic. For study of this system below 
the critical flutter speed, the motion must be assumed 
as damped. In this case the motion will be harmonic 
é“, modified by a typical damping factor e~™, 


This is equivalent to assigning a complex value to the 
flutter frequency : 


wo =o, + ir 


This is a legitimate procedure provided that the 
determinant elements themselves are accounted for 
under this assumption. In the flutter stability deter- 
minant each element is a combination of aerodynamic 
and mechanical terms. For the mechanical terms there 
is no error involved in assuming w complex. For the 
“noncirculatory’’ aerodynamic terms* there is no 
error in assuming w complex. However, in the “‘circu- 
latory’’ aerodynamic terms there is a function 


C(k) = F(1/k) + 1G(1/k) 


which is developed expressly on the assumption that 
the airfoil motion is harmonic (w is real). For such 
terms the assumption of a damped oscillatory motion 
is in error. However, the basic fact that the airfoil 
oscillates (even though damped) assures that the har- 
monic hypothesis upon which C(k) is calculated is ap- 
proximated fairly well. 

It might be noted that the approximation inherent 
in the aerodynamic effect of damped oscillations can 
actually be evaluated.*~® However, to retain sim- 
plicity the authors have not taken this effect into ac- 
count. Thus, there results a typical flutter analysis 
in which the structural damping required for the exist- 
ence of a flutter condition below the critical flutter 
speed are used with little or no modification. 

In the usual solution of the flutter polynomial the 
two unknowns chosen are g and w. The choice of g 
as one of the unknowns is arbitrary; instead, any one 
of several physical parameters (e.g., frequency ratio, 
aileron balance) could be taken as unknown. If now 
w = w; + 1d is allowed to be complex (i.e., the motion 
is assumed damped), the frequency w and the rate of 
decay \ may be chosen as the unknowns. Then the 
actual value g, of damping coefficient must be substi- 
tuted for g. The value of g, is detérmined from ground 
vibration test. The same roots Z of the characteristic 
flutter stability polynomial will then have the form 
shown previously with the new values for g and w: 
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Z=X+iY=(1+ = (1+ 


From this relation 


wo + ik = wV (1 + ige/(1 + ig) 


=> w[(1 + 1/2tga)/(1 + 


1 + (g?/4) 
w[1 + (i/2) (ga — 


since usually g, and g arelessthan0.1. Thus, 
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V= 347_MPH (TAS). 
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A = 


We 
Vx g) 

With these approximations there can be predicted the 
oscillatory frequency w, which should follow a mechani- 
cal disturbance of the airfoil at any speed up to the 
critical flutter speed and the decay exponent \ which is 
proportional to the logarithmic decrement of the 
damped oscillatory wave to be expected below the criti- 
cal flutter speed. 

The theoretical procedure is summarized as follows: 


(1) For a given 1/k, obtain the usual characteristic 
rootsZ = X +7Y. 

(2) Calculate g = Y/X. 

(3) Calculate v(m.p.h.) = w,(c.p.m.) (b,/14) [(1/k)/ 
VX]. 

(4) Calculate the flutter frequency, which is also the 
value taken for the frequency of damped oscillation: 
= = w. 

(5) Obtain g, from ground vibration test results. 

(6) Calculate \ = (wa/2WX) (g, — g). 

(7) Plot g vs. v(m.p.h.)/w.(c.p.m.), the standard 
g-V curve. 


AS 


(8) Plot A vs. v (m.p.h.). 

(9) Plot +e~™ as the envelope of a sinusoidal wave 
of frequency 

Then predicted decay characteristics of the motion of 
aircraft and wing below its critical flutter velocity can 
be compared with test results over a range of veloci- 
ties. 

The type of flutter analysis used by the authors in 
the subsequent example is that outlined in reference 
1—namely, a Rayleigh-type analysis based on un- 
coupled wing vibration modes. 


TEST 


The test procedure is as follows: 

(1) Obtain records of the oscillatory characteristics 
of a given aircraft wing after passage of a forcing im- 
pulse at a given true air speed (7) test- 

(2) Obtain the true oscillatory frequency w from test 
records. 

(3) Determine Bes from the logarithmic decrement 
of the test oscillations, as follows: 


(A)tese = (w/2m) loge + 1) 
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Test 2 __V= 443 

| 


| 


where A; is the maximum amplitude of the ith 
wave. 

(Test measurements should be made after the forcing 
impulse has ceased. In the results given in this paper 
a series of \ values for successive waves of each oscilla- 
tion were calculated and numerically averaged to ob- 
tain a “test value”’ of \.) 


(4) Plot (A)test VETSUS (2)test- 


Figs. 1-6 show the results of a series of tests made on 
a large four-engined aircraft of conventional design at 
the Republic Aviation Corporation. The flutter mode 
considered was unsymmetric, involving bending and 
torsion of the wings and roll of the entire aircraft. 
During level flight at various speeds the ailerons were 
abruptly deflected and released, causing unsymmetric 
oscillations of the wings. The test results were obtained 
on a recording oscillograph using electromagnetic 
accelerometer pickups located near the trailing edge 
of the wing in each outboard nacelle. The decay char- 


acteristics of the waves shown in Figs. 3-6 were ob- 
tained only for that portion of the motion which was 
free of the original forcing impulse of the ailerons. 


ij if 


The aileron rotation records (not shown) were found 
to damp out much more rapidly than those of the wing 
itself. Hence, theoretical prediction of the motion 
made on the basis of wing bending and torsion and air- 
plane roll (without aileron degree of freedom) was con- 
sidered adequate for comparison with test. 

The above outlined procedure for measuring decre- 
ments in flight is not claimed to be new. It has been 
used in various related manners both in this country 
and in Germany. It is possible that an oversimplified 
picture has been drawn without sufficient warning. It 
is realized that the modes of decay of an aircraft are 
complicated, even on the ground. They are, of course, 
dependent on the type of excitation or impulse given 
the aircraft. Experience indicates that for the case of 
a multiengined aircraft the unsymmetric flutter mode 
is usually the most critical one. Thus, in the particular 
example investigated, the excitation of the unsym- 
metric critical mode by a throw of the ailerons was 
feasible. Whether the scheme presented here can be 
generally used to predict flutter trends remains to be 
shown by other experimental results of considerable 


variety. 
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Stability of an Oscillating Airfoil in 
Supersonic Airflow’ 


M. V. BARTON 
The Unwersity of Texas 


SUMMARY 


The conditions for the determination of the stability of a two- 
dimensional airfoil oscillating harmonically with two degrees of 
freedom in a supersonic flow are developed. The stability is ex- 
pressed in terms of the position of the axis of rotation of the air- 
foil, the reduced frequency, the Mach Number, the ratio of the 
amplitudes of the torsional and flexural modes, and.the phase 
angle between the modes and evaluations are made for Mach 
Numbers ranging from 1.2 to 2.0 and for reduced frequencies 
from 0.1 to 2.0. The range of reduced frequencies and Mach 
Numbers for instability of a two degree of freedom system are 
indicated. It is shown that instability for single degree of free- 
dom torsional motion does not occur at a Mach Number greater 
than 1.58. The effects of the pivot position on the stability for 
various phase angles and amplitude ratios of the two modes of 
motion of the airfoil are discussed. — 

Tables of aerodynamic parameters for an oscillating airfoil in 
a supersonic flow are given as a supplement to previously pub- 
lished values. 


NOTATION 


free-stream velocity (ft. per sec.) 


M = Mach Number 

6 = Mach angle = sin~'(1/M) 

p = density of air (slugs per cu.ft.) 

a = angular displacement of airfoil about pivot (rad.) 

ao = amplitude of a 

y = transverse displacement of pivot (ft.) 

Yo = amplitude of y 

c = chord of airfoil (ft.) 

m = distance in fraction of chord aft from leading edge to 
pivot 

2 = circular frequency of oscillation of airfcil (rad. per 
sec.) 

k = reduced frequency = Qc/V cos? 6 

= amplitude ratio 

= phase angle 

Po = lift force per unit span (lbs. per ft.) 

Mo = aerodynamic moment per unit span about pivot 


a (Ibs. ft. per ft.) 
Pa, Ma*, M, etc. = aerodynamic parameters 


INTRODUCTION 


DETERMINATION of the air forces on an oscillat- 
ing airfoil in a supersonic flow is important because 
of the problem of flutter and aerodynamic instability 
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of high-speed aircraft and supersonic missiles. Several 
investigators including Possio,' Borbely,? Miles,* and 
Temple and Jahn‘ determined the air forces on an os- 
cillating airfoil in a two-dimensional supersonic flow 
by means of a linearized theory. This theory and the 
method of solution of the equations are fully discussed 
in the references cited and will not be considered here. 
Other investigators such as Garrick and Rubinow,”: ® 
Keller, Komm, and Johnson; and Barton and Poin- 
dexter,’ as well as some of those 1eferred to formerly, 
have made numerical evaluations ot the parameters 
involved. Some of these investigators, along with 
others, have applied the air-force calculations to the 
problem of flutter at supersonic speed in which the ef- 
fect of the structural parameters such as elastic axis 
position, structural damping, stiffness, and inertia 
have been considered. This report is concerned with 
the stability of a two degree of freedom system con- 
sisting of an airfoil oscillating torsionally and trans- 
versely simultaneously and does not consider the effect 
of the structural parameters on the motion. 

The motion of an airfoil oscillating harmonically 
with two degrees of freedom can be characterized by 
the amplitude ratio of the two components of motion 
and the phase angles between them. In an actual flut- 
ter situation the amplitude ratio and phase angle de- 
pend on the structural characteristics of the wing as 
well as the aerodynamic forces. This report does not 
consider the determination of the amplitude ratios and 
the amplitude ratios as such; it is only concerned with 
the stability of the system for specified values of these 
quantities. Since flutter cannot occur unless the motion 
is such that the air is supplying energy to the system, 
then, if the values of reduced frequency and Mach 
Number for this condition can be determined for speci- 
fied motions of the airfoil, the limits of the aerodynamic 
parameters in which flutter occurs can be established. 
In other words, the regions of the aerodynamic param- 
eters in which flutter can occur can be established in 
terms of the motion of the airfoil. 

The results obtained in this investigation are pri- 
marily useful for a better understanding of the phe- 
nomenon of supersonic flutter without the introduction 
of the complications of structural effects. It should be 
emphasized that the results and methods do not replace 
the usual type of flutter investigation in which the 
structural parameters are included, although the re- 
sults may act as a guide in selecting the range of values 
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Forces on airfoil. 


of Mach Numbers and reduced frequencies. to be con- 
sidered. 

The stability is investigated for Mach Numbers 
from 1.2 to 2.0 and for reduced frequencies of 0.1 to 
2.0. 


AERODYNAMIC FORCES 


The aerodynamic parameters used here have been 
evaluated by using the Borbely equations for the air 
forces on an oscillating airfoil in a supersonic flow. 
The assumptions and limitations of the theory have 
been thoroughly discussed in the references cited and 
will not be considered here except to repeat that the 
flow is assumed two-dimensional, an infinitely thin 
airfoil is used, boundary layer and friction effects are 
neglected, the airfoil performs small oscillations, and 
the results are not valid for Mach Numbers near unity. 

Consider an airfoil as shown in Fig. 1 oscillating har- 
monically in rotation about a pivot a distance yc behind 
the leading edge, while simultaneously the pivot is 
moving harmonically transversely to the undisturbed 
fluid flow direction. The aerodynamic lift and moment 
at the pivot per unit length of span are given by 


Po 


(P, + iP,*)a + (P, + iP,*) 
(1) 
My 


+iM,*)a+ (MM, +iM,*) 
Cc 


The values of the aerodynamic parameters, P,, and so 
on, depend on the pivot position, the Mach Number, 
and the reduced frequency. 

Since the airfoil is assumed to oscillate harmonically, 
the angular rotation and transverse displacement can 
be taken as 


a= 


y= (2) 


STABILITY CRITERION 


The work done by the lift and moment per cycle of 
oscillation is 


W= Podyt+ Moda 


where P, dy and M) da are scalar products. 
Substituting the values of Eqs. (1) and (2) into (3) 
gives _ 


W = Mat + P,) sin 6 + 
a 
yo \? * 
(M,* + P,*) cos¢] + P,*¢ (4) 


It is apparent that positive work is done and the sys- 
tem is unstable when 


M,* + = —P.) sin + (M,* + P,*)cos ¢] + 
0) 
2 
P,*>0 (5) 


The system is neutrally stable when the same quantity 
is zero and it is stable when this quantity is less than 
zero. 

For a single degree of freedom torsional motion, we 
let yo = 0 so that only rotation about the pivot is al- 
lowed.and the condition for torsional instability be- 
comes M,* > 0. 

For the case of single degree of freedom motion in 
flexure such that ay = 0, the condition for instability is 
P,* > 0. This mode of motion was found stable for 
all the values of Mach Number and reduced frequency 
investigated. 

It is sometimes convenient to be able to delineate the 
stable and unstable regions in terms of the Mach 
Number, reduced frequency, and pivot position, as 
well as the amplitude ratio and phase angle. This is 
accomplished by substituting the aerodynamic param- 
eters in terms of uw into the condition for neutral sta- 
bility. The substitutions to be made are: 


M,* = + + 


= Ms + + 

M,* = Ms + 

M, = Me + (6) 
= Pst+ 

P, = P, + uP, 

= P, 

Pp, = 


so that the condition for neutral stability becomes 


+ (Me — sin + (Ms + Ps) cos + 
ao 


=0 (7) 
aoe 


The values for the M/’s and P’s for various reduced 
frequencies and Mach Numbers are given in Tables 1 
through 5. 


These values are included here as a supple- 
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TABLE 1 


Mach Number 1.2 

k P2 Pi P, Ps Ms M M2 Ms Ms 
9.03056 0.1 3.4197 98.506 3,224.95 —62.624 —1,611.4 3,227.4 41.723 —13.41 2.279 49.209 
006111 0.2 —3.3992 —48.998 802.92 —31.046 —400.34 805.18 20.649 —6.633 2.263 24.414 
0.09167 O.3 3.3664 —32.385 354.43 —20.405 —176.11 356. 66 13.534 —4.342 2.236 16.065 
0.12222 0.4 —3.3205 —23.997 197.46 —15.003 —97.647 199.66 9.9091 —3.169 2.200 11.832 
0.15278 0.5 —3.2626 —18.905 124.86 — 11.696 —61.357 127.01 7.6837 —2.45 2.153 9.2476 
0.18333 O.6 —3.1927 —15.464 85.442 —9.4409 —41.675 87.540 6.1589 —1.951 2.098 7.4901 
0.21389 0.7 —3.1119 —12.969 61.716 —7.7895 —29. 833 63.748 5.0382 —1.581 2.0335 6.2090 
0.24444 0.8 —3.0211 —11.071 46.351 —6.5171 —22.176 48.311 4.1714 —1.289 1.9614 5.2279 
0.27500 O.9 —2.9214 —9.5723 35.847 —5.5007 — 16.954 37.731 3.4764 —1.050 1.8822 4.4503 
0.30556 1.0 —2.8131 —8.3576 28.369 —4.6668 —13.247 30. 166 2.9047 —0.8491 1.7968 3.8147 
0.33611 1.1 —2.6978 —7.3529 22.869 —3.9676 — 10.530 24.575 2.4238 —0.6745 1.7061 3.2937 
0.36667 1.2 —2.5769 —6.5084 18.716 —3.3718 —8.4905 20.327 2.0141 —0.5189 1.6110 2.8527 
0.39722 1.3 —2.4511 —5.7902 15.515 —2.8589 ~—6.9278 17.027 1.6604 —0.3802 1.5125 2.4787 
0.42778 1.4 —2.3216 —5.1739 13.004 —2.4121 —5.7115 14.416 1.3544 —0.2527 1.4119 2.1591 
0.45833 1.5 —2.1896 —4,6411 11.006 —2.0210 —4.7530 12.315 1.0844 —0.1360 1.3093 1.8851 
0.48889 1.6 —2.0564 —4.1786 9.3969 — 1.6770 —3.9893 10.604 0.84878 —0.0271 1.2066 1.6497 
0.51944 1.7 —1.9229 —3.7761 8.0877 —1.3731 —3.3760 9.1921 0.64191 0.0745 1.1043 1.4475 
0.50000 1.8 —3.4246 7.0134 —1.105 —2.8806 8.0168 0.46016 0.1698 1.0034 1.2743 
0.58056 1.9 —3.1176 6.1252 —0.8668 —2.4781 7.0300 0.300813 0.2596 0.90468 1.1264 
0.61111 2.0 —2.8497 5.3868 —0.6573 —2.1499 6.1957 0.161441 0.3437 0.80898 1.0009 

TABLE 2 
Mach Number 1.4 

k Ps P; M; My M: Ms Ms 
0.04898 0.1 —1.061 —41.622 850.12 —0.841 —424.90 850.83 0.557 19.957 0.70717 20.796 
0.09796 0.2 —1.055 —20.730 211.98 —0.380 —105.81 212.69 0.2454 9.959 0.70236 10.339 
0.14694 0.3 —1.0453 —13.734 93.818 —0.208 —46.729 94.512 0.128 6.620 0.69446 6.8276 
0.19592 0.4 —1.0316 —10.210 52.463 —0. 1086 — 26.054 53.147 0.0582 4.9444 0.68378 5.0530 
0.244909 0.5 —1.0142 —8.0766 33.325 —0.0392 — 16.487 33.994 0.0082 3.9352 0.66964 3.9744 
0.29388 0.6 —0.99323 —6.6398 22.933 0.016 —11.294 23.586 —0.0314 3.2600 0.65298 3.2445 
0.34286 0.7 —0.96896 —5.6019 16.674 0.0610 —8. 1677 17.308 —0.0649 2.7758 0.63363 2.7146 
0.39184 O.8 —0.94161 —4.8144 12.616 0.1002 —6. 1426 13.228 —0.0941 2.4111 6.61187 2.3108 
0.44082 0.9 —0.91135 —4.1947 9.8394 0.135 —4.7586 10.427 —0.1202 2.1271 0.58791 1.9919 
0.48980 1.0 —0.87855 —3.6936 7.8578 0. 167 —3.7727 8.4197 —0.1439 1.900 0.56194 1.73: 
0.53878 1.1 —0.84344 —3.2798 6.3969 0.195 —3.0473 6.9311 —0.1653 1.714 0.53422 et 
0.58776 1.2 —0.80642 —2.9324 5.2904 0.2217 —2.4994 5.7954 —0.1849 1.560 0.50505 1.3388 
0.63673 1.3 —0.76767 336 4.4337 0.2454 —2.0771 4.9085 —0. 2027 1.431 0.47466 1.1857 
0.68571 1.4 —0.72761 25 3.7587 0.2670 —1.7455 4.2022 —0.21880 1.321 0.44337 1.0545 
0.73469 1.5 —0.68654 —2.1624 3.2186 0.2864 —1.4814 3.6301 —0.23335 1.228 0.41140 0.94136 
0.78367 1.6 —0.64479 —1.9706 2.7803 0.30432 — 1.2687 3.1594 —0.24631 1.148 0.37909 0.84340 
0.83265 1.7  —0.60277 25 2.4209 0.32012 —1.0954 2.7677 —0.25778 1.079 0.34666 0.7584: 
0.88163 1.8 —0.56066 2.1235 0.33435 —0.95318 2.4378 ~-0. 26795 1.0191 0.31441 0 
0.93061 1.9 —0.51881 1.8751 0.34674 —0.83548 2.1575 —0.27656 0.96755 0.28256 
0.97959 2.0 —0.47756 “S104 1.6661 0.35752 —0.73751 1.9175 —0 38387 0.92301 0.25141 

TABLE 3 
Mach Number 1.6 

k P2 P; Ps Ps; Ms Ms Mi M2 Ms Ms 
0.06044 0.1 —26.469 438 13 4.5848 —219.01 438.48 —3.0569 17.813 0.34865 13.228 
0.12088 0.2 —13.196 2. 3060 —54.608 109.69 —1.5401 8.8904 0.34636 6.5849 

0.3 51562 —8.7541 1.5527 —24.164 48.798 — 1.0388 5.9104 0.34258 4.3578 
0.4 —0.50905 —6.5211 1.1802 —13.511 27.485 —0.79170 4.4153 0.33735 3.2351 
0.5 —0.50071 —5.1708 0.96040 —8.5807 17.616 —0.64633 3.5149 0.33072 2.55 

36263 0.6 —0.49067 —4.2652 0.81671 —5.9040 12.253 —0.55149 2.9121 3 2 
0 12306 0.7 —0.47904 —$.6137 0.71613 —4.2913 9.0163 —0.48567 2.4793 3 1 
0.48350 0.8 —0.46587 0 64249 —3.2462 6.9127 —0.43751 2.1534 0.30297 1. 
0.54394 0.9 —0.45131 —2.7286 58694 —2.5309 5.4678 —0.40149 1.8990 0.29140 1.3121 
0.60438 1.0 —0.43546 —2.4152 54339 ~2.0207 4.4316 —0.37344 1.6945 0. 27887 1.15 
0.66481 1.1 —0.41847 — 2.1565 0: 50910 —1.6444 3.6623 —0.35129 1.5272 0.26542 
0.72525 1.2 —0.40042 — 1.9397 2.8237 0.48114 596 3.0748 —0.33347 1.3875 0.25121 
0.78569  —0.38156 —1.7552 2.3791 0.45821 2.6156 —0.31887 1.2696 0.23638 
0.84613 1.4 —0.36197 — 1.5966 2.0279 0.43896 2.2489 —0. 30666 1. 1687 0.22104 
0.90656 1.5 —0.34180 — 1.4590 1.7459 0.42271 1.9512 —0.29635 1.0820 0.20530 
0.96700 1.6 —0.32122 — 1.3388 1.5164 0.40873 1.7057 ~0. 28750 1.0067 0. 18932 
1.02744 1.7 —0.30037 — 1.2333 1.3276 0.39653 1.5007 —0. 27973 0.94099 0.17318 
1.08788 1.8 —0.27939 —1.1403 1.1705 0.38574 1.3276 —0. 27283 0. 88349 0.15706 
1.148381 1.9 —0.25845 — 1.0580 1.0388 0.37610 1.1799 —0. 26657 0.83296 0.14104 0.45687 
1:20875 2.0 —0.23766 —0.98504 0.92742 0.36734 1.0527 —0. 26084 0.78846 0.12525 0.42113 

TABLE 
Mach Number 1.8 

k P2 Pi Ps Ps Ms Ms M M2 = Ms 
0.6914 0.1 --0.29773 —19.314 279.47 5. 353 —139.71 279.66 —3.5693 15.007 0.19841 9.6530 
0.13827 0.2 —0.29615 —9.6345 69.778 2. 6829 ~34.860 69.975 — 1.7897 7.493 0.19714 4.8099 
0.20741 0.3 —0.29346 —6 3986 30.949 1.7952 —15.446 31.144 —1.199 4.9835 0.19502 3.1881 
0.27654 0.4 —0.28983 —4.7735 17.359 1.3534 —8.6511 17.551 —0.9043 3.7256 0.19211 2.3722 
0.34568 0.5 —3.7931 11.069 1.090 —5.5066 11.258 —0.7292 2.9681 0. 18834 1.8784 
0.41481 0.6 —3.1351 7.6537 0.9151 —3.7991 7.8375 —0.6131 2.4613 0. 18385 1.5462 
0 18395 0.7 —2.6619 5.5947 0.7915 —2.7701 5.7732 —0.5312 2.0981 0.17860 1.3066 

: 0.8 —2.3043 * 4.2591 0.6994 —2.1028 4.4317 —0.4704 1.8243 0.17268 1.1248 

5 0.9 257 —2.0241 3.3441 0.6287 —1.6459 3.5103 —0.42351 1.6106 0.16614 0.98205 
0 60136 1.0 —0.24826 — 1.7983 2.6904 0.5725 —1.3197 2.8494 — 0.38648 1.4392 0. 15902 0.86667 
0.76049 1.1 —0.23858 —1.6122 2.2073 0.5269 ~1.0789 2.3587 —0.35654 1.298 0.15140 0.77147 
0.82963 1.2 —0.22840 —1.4564 1.8405 0.4896 —0.89624 1.9839 —0.33198 1.181 0.14333 0.69160 
0.89877. 1.3 —0.21759 ~ 1.3238 1.5557 0.45799 —0.75465 1.6906 —0.31124 1.082 0.13485 0.62366 
0.96790 1.4 —0.20644 —1.2100 1.3304 0.43142 —0. 64283 1.4566 —0. 29379 0.9969 0.12609 0.56536 
1.08704 1.5 —0.19486 —1.1111 1.1492 0.40839 —0.55311 1.2663 —0.27870 0.9232 0.11705 0.51485 
1.10617 1.6 —0.18303 — 1.0247 1.0015 0.38841 —0.48016 1.1094 —0. 26559 0.85927 0.10784 0.47086 
1.17531 1.7. —0.17102 —0. 94873 0.87971 0.37088 —0.42016 0.97823 —0.25410 0.80322 0.09851 0.43235 
1.24444 1.8 —0.15889 —0.88157 0.77814 0.35534 —0.37030 0.86731 —0.24387 0.75391 0.08917 0.39858 
1.31358 1.9 —0.14675 —0.82198 0.69267 0.34144 ~—0.32853 0.77252 —0.23473 0.71031 0.07984 0.36887 
1.38272 2.0 0.13463 — 0.76896 0.62015 0.32894 —0, 29324 0.69076 —0. 22645 0.67166 0.07061 0.34272 
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TABLE 5 


Mach Number 2.0 


k Ps Pi Ps Ps Ms Ms M: Ms Ms 
0.07500 0.1 0.19210 15.386 205.22 5.1330 102.59 205.35 —3.4226 12.824 0.12802 7.6905 
0.15000 0.2 —0.1910¥ —7.6788 51.256 2.5705 —25.612 51.383 —1.7143 6.405 0.12720 3.8345 
0.22500 0.3 —0.18939 5.1033 22.745 1.7169 —11:357 22/871 —1.1455 4.2615 0.12585 2.5445 
0.30000 0.4 —0.18706  -3.8111 12.766 1.2914 —6.3675 12.891 -0. 3.1875 0.12398 1.8961 
0.37500 0.5 —0.18408  —3.0322 8.1486 1.0369 —4 8.2702 —0.69277 2.5414 0.12159 1.5045 
0.45000 0.6 —0.18047 2.5103 5.6402 0. 86786 —2.8048 5.7590  —0.58023 2.1093 0.11873 1.2414 
0.52500 0.7 —0.17629 —2.1353 4.1281 0.74769 — 2/0489 4.2434 —0.500 1.7996 0.11540 
0.60000 0.8 —0.17154  —1.8523 3.1470 0.65792 — 1.5587 3.2587 —0.44085 1.5664 0.11161 0.90847 
0.67500 0.9 —0.16625 —1.6308 2.4748 0.58856 — 112228 2.5821 —0.39473 1.3845. 0.107420. 79595 
0.72 1:0 0.16050 1.4525 1.9943 0.53326 —0:98294 2.0971 —0.35808 1.2384 0.10285 0.70516 
0.82500 1.1 0.15431 —1.3057 1.6390 48830 —0. 805 1.7369  —0.32835 1.1187 0.097953 0.63037 
0.9 1.2 -—0.14773  —1.1829 1.3694 0.45106 —0.67124 1.4621 —0.30370 1.0188 0.092769 0.56767 
0.97570 1.3 —0.14079 —1.0785 1.1597 0.41971 0.56682 1.2470  —0.28297 0.93404 0.087307 0.51437 
1.05000 1.4  —0.13360  —0.98886 0.99370 0.39300 0.48428 1.0754 —0.26533 0.86158 0.08166 0.46861 
1.12500 1.5 0.12615 —0.91104 0.86012 0.36993 0.41792 0.93594 —0.25010 0.79887 0.075818 0.42894 
1.20000 1.6 —0.11848  —0.84294 0.751 0.34985 —0.36388 0.82093 —0.23683 0.74420 0.069848 0.39435 
1.27500 1.7 0.11069  —0.78300 0.66102 0.33217 —0.31933 0.72482 —0.22514 0.69619 0.36402 
1.35000 1.8 —0.10281 —0.72997 0.585 0.31648 —0:28223 0.64354. —0.21479 0.65381 «0.057712 «0.33733 
1.42500 1.9 —0.09489  —0.68281 0.52244 0.30249 —0.25103 0.57406 —0.20550 0.61623 0.051615 0.31374 
1.50000 2.0  —0.08700  —0.64077 0.46861 0.28988 —0.22462 0.51419 -0.19718 0.58278 0.045576 0.29290 
ome i but merely moves the range of pivot positions forward 
spot so that for an amplitude ratio of unity, the unstable 
Pe mony region extends for a range of pivot positions of about 
Tr \ one-half to a full chord ahead of the wing. It should 
we \) also be noticed that the maximum value of the reduced 
oz UNSTABLE frequency for which instability can occur does not 
od f, change with the amplitude ratio for this particular 
-g4+0 a _ ra phase angle and Mach Number. If the phase angle 
« “aoe oss a Gaal = were 180° instead of zero, then a similar set of curves 
z —T would result except that they would be displaced aft 
e of the leading edge of the wing as the amplitude ratio is 
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° increased. 
pa A similar set of curves for a phase angle of 270° and 
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ment to previously published values for somewhat dif- | 
ferent Mach Numbers and reduced frequencies. > 
a 
The lines of neutral stability for a Mach Number of 1.4 - / 
and for various values of reduced frequency are shown sal STABLE 
in Figs. 2 and 3. Fig. 2 shows the stable and unstable V4 
regions in terms of the reduced frequency and pivot P L 
position for various values of the amplitude ratio and / 
for zero phase angle. The dashed-line curve is for / 
yo/acc = 0, and it therefore corresponds to the single i / 
degree of freedom torsion case. It is evident that tor- | 
sional instability can occur at this Mach Number for | 
pivot positions ranging from slightly ahead of the lead- a 


ing edge of the wing to slightly behind the mid-chord. 
Increasing the magnitude of the amplitude ratio does 
not change the size or the shape of the unstable region 
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Fic. 3. Stable and unstable regions for various values of ampli- 
tude ratio at Mach Number 1.4 and 270° phase angle. 
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STABILITY OF AN OSCILLATING AIRFOIL 


the curves separating the stable and unstable regions 
change size and shape with changing amplitude ratio. 
The unstable regions cover a large range of pivot posi- 
tions for small values of the reduced frequency, and in 
general this range increases with an increase of ampli- 
tude ratio. The maximum value of the reduced fre- 
quency for which instability can occur is greater for 
every nonzero value of yo/aoc investigated than for the 
case of torsional oscillations for which yo/aoc = 0. 

The regions of instability for a phase angle of 90° and 
a Mach Number of 1.4 are of little interest. In fact, 
instability only occurs for extremely small values of the 
amplitude ratio and then the range of pivot positions is 
well within the limits determined for the other phase 
angles previously discussed. 


LIMITING VALUES OF REDUCED FREQUENCY FOR 
INSTABILITY 


The curves of Figs. 2 and 3 indicate that for a given 
Mach Number and amplitude ratio there are limits to 
the values the reduced frequency can have for insta- 
bility regardless of the pivot position. These limiting 
values of the reduced frequency for a specified airfoil 
motion and Mach Number can be determined by mak- 
ing the discriminant of Eq. (7) equal to zero. Since the 
condition for neutral stability is a quadratic in y, only 
real roots are obtained when the discriminant is posi- 
tive. The values of the reduced frequency required to 
make the discriminant zero corresponds to the limiting 
high and low values of the reduced frequency for neu- 
tral stability. 

The value of the discriminant, D, of Eq. (7) is 


D = My? — 4M,P, + — — 
a 


P;) cos ¢ — (Ms — Ps) sin ¢] — sin? @ 


ac 
It can be shown that 
M, — Ms —- P; = 0 
so that 


~ + — P,) sind — 
a 
2 
1(2*) Py (8) 


The values of the reduced frequencies necessary to 
make the discriminant zero for various values of ampli- 
tude ratios, phase angles, and Mach Numbers are 
shown in Figs. 4 and 5. 

In Fig. 4 the values of k for neutral stability are 
plotted radially and the values of the phase angles are 
plotted angularly for various values of the amplitude 
ratio and for a Mach Number of 1.4. The dashed-line 
circle corresponds to the maximum value of the reduced 
frequency for which single degree of freedom torsional 


375 


EXTRAPOLATED VALUES 


Fic. 4. Stable and unstable regions for various amplitude 
ratios at Mach Number 1.4. 


EXTRAPOLATED VALUES 


Fic. 5. Stable and unstable regions for various Mach Numbers 
and for an amplitude ratio of 0.1. 


instability can occur. It is interesting to note that for 
an extremely small transverse motion such that yo/aoc = 
0.01, the maximum value of k corresponding to in- 
stability does not differ much from the torsional case. 
However, when there is a slight transverse motion 
combined with.the rotation, there is a minimum value 
of k below which instability cannot occur. These lower 
values of k for neutral stability are indicated by the 
small loop for values of phase angle in the neighborhood 
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of 90°. As the amplitude ratio increases in value, the 
maximum value of the reduced frequency for neutral 
stability increases and the maximum occurs at different 
phase angles. Thus, for yo/acc = 0.1, the maximum 
value of k corresponds to a phase angle of 270°, whereas, 
for yo/acc = 1.0, the maximum values seem to be 
around 205° and 335°. 

The effect of variation of Mach Number on the sta- 
bility is shown in Fig. 5. For an amplitude ratio of 
0.1 it is evident that increasing the Mach Number de- 
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creases the values of the reduced frequency for which 
instability can occur. q 

The values of the reduced frequency and Mach Numa 
ber for unstable flow in the special case of only tom 
sional oscillations is shown in Fig. 6. It should be 
noted that for this single degree of freedom motion, ins 
stability cannot occur at a Mach Number of over 1.58 
As the Mach Number is decreased from this value, thé 
range of reduced frequencies for instability of the sys: 
tem increases. 
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